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During  the  fall  semester  of  the  academic  year  1952-1953 
Dr.  C.  J.  Bouwkamp  of  the  Philips  Research  Laboratories,  Eindhoven, 
Netherlands  served  as  a  member  of  the  Washington  Square  College  Mathe- 
matics Research  Group  of  New  York  University,  Dr.  Bouwkamp  is,  of 
course,  well  known  for  his  research  in  electromagnetic  theory  and  re- 
lated fields.  In  recent  months  Dr.  Bouwkamp  has  been  surveying  and 
assessing  the  literature  on  the  subject  of  diffraction.  The  present 
report  concentrates  on  papers  related  to  the  subject  of  diffraction  by 
apertures.  However,  the  report  does  not  merely  relate  the  efforts  of 
numerous  authors,  but,  where  necessary,  criticizes,  emends,  and  pixrsues 
the  consequences  of  their  work. 

This  study  of  the  problem  of  diffraction  illus- 
trates very  clearly  the  need  for  mathematical  accuracy  in  the  formulation 
and  solution  of  physical  problems.  Because  the  diffraction  problem  is  a 
difficult  one,  many  authors  have  made  plausible  physical  ass'jimptions  which 
simplify  the  procedure.  Others  have  introduced  reasonable  mathematical 
approximations  at  some  stage  of  their  work.  However,  further  mathematical 
analysis  has  shown  many  of  these  solutions  to  be  inadequate;  some  are  in 
error  for  certain  ranges  of  wavelengths  or  in  certain  regions  of  space, 
and  some  lead  to  mathematical  ambiguities  or  even  inconsistencies,  Dr, 
Bouwkamp 's  extensive  investigations  of  these  papers  and  his  passion  for 
accuracy  (which  is,  however,  tempered  by  a  full  appreciation  of  the  pre- 
sent state  of  applied  mathematics)  have  made  apparent  the  limitations  of 
many  of  the  solutions  that  have  been  offered.  His  work  exemplifies  a 
valuable  principle  of  scientific  research,  namely,  that  a  published  paper 
tells  us  where  to  concentrate  our  doubts. 

Those  of  us  who  were  associated  with  Dr.  Bouwkamp  during  his 
stay  at  New  York  University  were  rewarded  both  personally  and  scientifically. 


Morris  Kline 
Project  Director 
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ABSTRACT 

A  number  of  recent  developments  in  the  theory  of  diffraction 
of  electromagnetic  waves,   particularly  those  dealing  with  apertu37es 
in  plane  conducting  screens,   are  reviewed.     The  subjects  treated 
include  modifications  of  Kirchhof f ' s  theory,  the  theory  of  small 
apertures,  Babinet's  principle  for  plane  obstacles,  variational 
principles,   and  singularities  at  sharp  edges. 

For  completeness,  a  discussion     from  an  alternative  view- 
point of  the  problem  of  diffraction  by  an  aperture  by  Professor 
N,  Marcuvitz  has  been  included  in  this  report. 
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I.   Introdnetion 

The  theory  of  diffraction  has  three  major  fields  of  application:  (l)  optica, 
(2)  radio-wave  propagation,  and  (3)  acoTistics.   In  (2)  and  (3)  the  wavelengths  con- 
sidered are  us-ually  of  the  same  order  of  magnitude  ae  the  diffracting  obstacle, 
while  in  the  case  (l)  the  wavelengths  are  tisually  small  compared  to  the  o'bstacle. 
A  further  difference  "between  these  three  fields  is  that  (l)  and  (2)  involve  es- 
sentially vectorial  problems,  while  the  problems  involved  in  (3)  are  mainly  scalar. 
However,  in  most  applications  of  diffraction  theory  to  classical  optics,  light 
is  considered  as  a  scalar  wave  phenomenon  (polarization  effects  are  ignored). 
For  example,  calculations  on  diffraction  by  apertures  are  ustjally  based  on 
Kirchhoff's  mathematical  formulation  of  Huygens •  principle.  Experiments  have 
shown  that  this  is  JTistifiable  when  the  wavelength  is  small  in  comparison  to  the 
size  of  the  aperture.  Polarization  cannot  be  ignored  in  radio-wave  propagation, 
where  the  wavelength  is  of  the  same  order  of  magnitude  as  the  apert-'ore.   ftie  way 
to  cover  this  is  by  using  an  electromagnetic  equivalent  of  the  scalar  Huygens- 
Kirchhoff  form^ila.  This  scalar  formula  may  be  applied  to  any  of  the  six  rec- 
tengiilar  components  of  the  electric  and  ma,°Tietic  vectors.   In  order  that  the 
six  wave  functions  so  obtained  should  satisfy  Maxwell's  equations  we  have  to 
introduce  certain  contour  integrals  along  the  rim  of  the  aperture  (Kottler). 

The  theory  of  Kirchhoff  and  Kottler  are  poor  substitutes  for  rigorous 
diffraction  theory  (wave  equation  plus  boundary  conditions)  in  the  qijasi-optical 
ranre  because  they  do  not  correctly  describe  the  field  in  the  vicinity  of  the 
aperture  and  the  edge.   In  the  extreme  case  of  very  long  waves  they  entirely 
fail  to  predict  the  correct  order  of  magnitude  of  the  field  far  from  the  aperture 
(Rayleigh). 

The  purpose  of  this  report  Is  to  comment  on  some  of  the  new  developments 
in  diffraction  theory.  Various  modifications  of  the  Kirchhoff  theory  have  re- 
cently been  proposed.   Bayleigh's  potential-approach  has  been  extended  to 
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hig'her-order  approxltnatlons.  The  Integral-eq-oation  technique  has  heen  developed 
extensively,  and  variational  methods  have  shown  their  usefulness  in  a  great  number 
of  prohlems.  Also,  the  rigorous  form  of  Babinet's  principle  in  plane  obstacle 
diffraction  theory  has  been  obtained.  New  insight  into  the  character  of  singu- 
larities at  sharp  edges  has  profoundly  influenced  many  aspects  of  diffraction 
theory. 

A  number  of  important  recent  developments  have  not  been  treated  in  these 
lectures.  Among  these  we  can  mention  the  exact  solutions  recently  obtained  for 
diffraction  by  circular  spertiires  and  disks,  and  the  Wiener-Hopf  technique  which 
has  proved  its  power  in  the  solution  of  certain  waveguide  problems. 

Only  steady-stste  problems  will  be  discussed.  The  time  factor  is  under- 
stood to  be  exTi (- itet ) .  For  a  general  introduction  into  diffraction  theory, 
which  includes  descriptions  of  the  early  work  by  Kirchhoff ,  Kottler  and 
Eayleigh,  see: 

Baker  and  Copson,  The  Mathematical  Theory  of  Huygens '  Principle.  Oxford, 

Clarendon  Press,  1950. 

Sommerfeld,  Vprlesungen  uber  theoretische  Physik.  vol.  k,   Optlk  , 
Wiesbaden,  Dieterich  Verlag,  1950. 


II.  Klrchhoff 3  Theory  of  Diffraction 

Let  ^  "be  a  screen  of  vanishing  thickness  covering  a  finite  part  of  the 

plane  «  =  0.   Consider  a  system  of  sources  in  the  left  half-space  (z  <  O),   If 

the  screen  were  absent,  these  sotirees  would  produce  a  wave  field  u  (P)  at  the 

point  P.  The  actual  field  u(P),  produced  when  the  diffracting  screen  is  present. 

Is  the  sum  of  u  (P)  and  u. (P),  where  u, (P)  is  the  diffracted  field  due  to  the 
0        0.  a 

secondary  sources   on  ^__.     By  Green's   theorem 


r   \        ,        ikr  Ikr    :,      1 


where  the  integration  is  over  hoth  faces  of  2_»  ^^^   r  1-s  the  distance  from  P  to 
d  /  ,  "P-urther,  h/hv  denotes  differentlption  with  respect  to  the  integration- 
point  coordinates  in  the  direction  of  the  normal  to  2_  drawn  into  free  space. 
Klrchhoff  made  the  following  assiunptions : 

du  ^\ 
(i)   u  =  u  ,  r—  =  r-—   on  S  (illuminated  side  of  screen) 
0  ov       ov 

(ii)  u  =  0  ,  ^  =  0         (on  dark  side  of  screen). 


The  total  field  then 'becomea ,   In  Klrchhoff  s  approximation, 

dZ. 


,  r  1   ikr  du        i    ikr 

(2.1)    i.(p) . u„(p) .y^^'—  ^  ..^1- (^) 


where  now  the  integration  Is  only  over  the  illuminated  part  S  of  X.  i  wliile  n 

refers  to  the  normal  of  S  drawn  Into  the  shadow  region  (z  >  O). 

Serious  ohjections  can  he  raised  against  Klrchhoff 's  theory*-  -^  .   In  fact,  as 

we  let  P  approach  a  point  Q  on  the  screen  /  ,  equation  (l)  falls  to  reproduce  the 

asstimed  values  u  and  9u  /hn.     This  can  "be  shown  with  the  use  of  the  following 
0      o 


theorem : 


^   /"    g    ikr  1  a  /   e^^         ^  J  P  ->  Q*" 


-  il  _ 


Then  the  limiting  values  of  K  (P)  are 


.i-h. 


ikr  ^n 


f^-^'  .^*e;^V  sr  '21, 


1    ? 
where  C  =  —  or  -r  according  as  P  is  on  the  darlc  or  on  the  illuminated  side  of 


/  .  Consequently,  only  if 


-      ,   /  ITcr  au 


will  the  limiting  values  of  K  he  identical  with  the  assumed  values  u  .  However, 

8  0 

this  condition  cannot  he  fulfilled  for  arbitrary  S  and  u  ,  as  can  easily  he  seen 
if  we  take  two  screens,  one  inside  the  other,  and  subtract  their  field  eqTiations; 
Then,  using  the  above  condition,  we  obtain 

ikr  du 


/   ikr  du 


where  §   (the  ••difference"  between  the   two  original   screens)  and  u    are  arbitrary. 

au 

It  would  then  follow  that  r— ^  =  0  on  S,  But  u  was  arbitrarily  chosen  and  there- 
^  3n  0 

ou 

fore  T— ^  would  not  necessarily  be  eaual  to  zero.  Hence  we  have  a  contradiction, 
on 

This  shows  that  Kirchhoff s  procedure  is  not  self-consistent.  The  same  conclusion 

follows  from  a  consideration  of  the  limiting  values  of  9K  /Bn,  which  are 

s 

9n      ,   -2    ;v2    .    /   ^ikr 

Trenerally  it  can  be  said  thR,t  the  reason  for  the  inconsistency  of  Kirchhoff  s 
theory  is  that  u  and  3u/au  cannot  simultaneously  be  prescribed  on  2Z  since  the 
equation  /^u+ku=  0  is  elliptic. 

If  Kirchhoff 's  boundary  conditions  on  ^  were  exact,  then  u  and  Bu/3n 
would  jump  by  the  amoun*^.s  u  and  9u  /Bn  respectively  across  S^   In  fact,  these 
jijmps  are  produced  by  K  ,  as  may  at  once  be  verified  from  expressions  (2)  and 


(3)  for  the  limiting  values  on  S.  This  is  in  accordance  with  Kottler^s  inter- 
he  : 
[2] 


pretation  of  Kirchhoff 's  formula  (l)  in  that  K  is  the  rigorous  solution,  not 

8 


of  a  boundary- value  problem,  but  of  a  saltus  problem' 
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We  shall  now  discuss  the  Complementary  problem  of  diffraction  "by  an  infinite 

plane  wcreen  with  a  finite  aperture  A,  Assuming  the  same  primary  field  u  ,  we 

o 

apply  equation  (l)  where  now  S  means  the  complement  of  A,  To  avoid  the  slight 

difficulty  thflt  S   is  no  longer  finite,    it  may  "be  necessary  to  assume  that  the 

Impglnpry  pert   of  k  is  positive.     Now  the  integral  over  S  is  equal  to  the  integral 

over  S  +  A  minus  the  integral  over  A,     The  Integral  over  S  +  A  equals  -u   (P) 

o 

if  2  >  0  and  equals  zero  if  z  <  0.  Hence  the  total  field  behind  the  aperture,  In 
Trirchho-*'f 's  sense.  Is 


a 


,   /•  f  ikr  du       :.    ikr  1 


(2.5)    ..(P)=  ujp).  ^ /{l^  ^  .  „^  I-  (1_)  J  a  Z  . 


and  the  total  field  in  front  of  the  aperture  is 

iter  9u      ^    iter 

.<' 

where  n  is  the  normal  of  A  drawn  into  the  shadow  region  (z  >  O). 

The  analytic  continuation  of  the  function  (h)   through  the  aperture  is  precisely 
the  function  (5).  and  vice  versa.  There  are  no  discontinuities  in  the  aperture, 
where  we  have 


.         .   f    iter  5u 


(2.6) 

aK^(P)    ,  3u^(P)      o  ,  d^    ,2      /   ,^^ 

Equation  {^)   is  usually  derived  in  a  different  manner,  namely  by  applying 
Green's  theorem  for  the  half-space  z  >  0  and  assuming  that  u  and  du/3n  are  zero 
on  the  dark  face  of  the  screen  and  that  in  the  aperture  they  are  equal  to  the 
unperturbed  values.  However,  Kirchhoff 's  original  method  is  preferable  since 
it  fivoids  the  difficulty  that  (ii)  does  not  reproduce  the  values  assumed  in  the 
aperture  but  rather  the  values  (6). 

Tor  complementary  problems  (A  =  S),  It  follows  from  (l),  (^)  and  (5) 
thPt 
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^o  +  ^»  =  ^^  f  or  z  >  0 

a    8    0 

K  *  K  -  ^  for   z  <  0. 

as     0 

This  is  Bflhinet's  principle  in  the  sense  of  Kirchhoff  s  theorj'^  •', 

It  has  often  been  suggested'-  ■■  that  Kirchhoff 's  formule  gives  the  first 
term  of  an  acc-urate  solution  of  a  ho-undary-value  problem  by  successive  approxi- 
mations. This  was  disproved  by  Franz"-  ■'  and  Schelkunoff  "■  •' .   In  addition, 
Franz  derived  equations  (^)  and  (5)  as  follows.   (A  similar  Interpretation  was 

given  by  Schelkunoff.)  A  zero-order  approximation  is  u  =  u  (z  <  O),  u  =  0 

o 

(z  >  O).  This  agrees  with  the  boundary  condition  at  the  black  screen,  but  the 
wave  equation  is  violated  in  the  aperture.  The  wave  equation  is  then  restored 
by  a  correction  term  arising  from  the  secondary  sources  in  the  aperture.  This 
term  equals  the  right-hand  side  of  equr-tion  {h)   for  all  P  (z  <  O).  By  elaborating 
the  new  interpretation,  Fre-nz""  ■"  devised  an  interesting  diffraction  theory  which 
is  applicable  for  all  wavelengths.   In  the  author's  opinion,  however,  it  is  still 
an  open  question  whether  Franz's  theory  will  eventually  become  of  value  in  practice, 


III.  Modified  Kirehhof^  Theory 

Modifications  of  Kirchhoff 's  theory  have  recently  been  proposed  for  planar 
diffraction  problems*.  One  aim  of  this  modified  Kirchhoff 's  theory  was  to  provide 
a  means  for  distinguishing  between  the  two  principal  boundary-value  problems: 
(l)  scalnr  wave  function  vanishing  on  a  soft  screen  and  (II)  normal  derivative  of 
scalsr  wave  function  vanishing  on  a  rigid  screen  (in  the  terminology  of  acousfecs). 
The  modified  theory  makes  use  of  the  two  principal  Green's  f-unctions  of  the  half- 
space,  which  are  known  explicitly. 

*  See, for  example,  Bouwkamp,  Thesis,  Groni''^en,  19^1. 
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let  ■ai'P)   denote  a  wave  function  that  is  regular  for  z  >  0.  Then,  in  the 
half-space  z  >  0,  we  have  Eayleigh's  formulea: 

(3.1)        1,:  u(P).-i/2^  U  dE. 

,    r     ^        Ikr 

(3.2)         V  u(p)=  i/u^c^)''!:. 

where  n  points  into  the  half-space  z  >  0,  and  the  integration  is  over  z  =  +0. 

EquPtlons  (l)  and  (2)  follow  at  once  from  an  application  of  Green's  theorem  using 

the  two  principal  Green's  functions  of  the  wave  equation  for  the  half-space. 

In  applying  Heyleigh's  formulas  to  diffraction  problems  we  shall  assume,  as 

in  the  usual  version  of  Kirchhoff 's  theory,  that  u  and  Su/9n  may  "be  replaced  "by 

the  corresponding  e-eometricsl-optics  values:      Immedietely  "behind  the  screen  u  and 

3u/^n  pre  tsken  to  "be  zero,  and  in  the  aperture  they  are  replaced  "by  u  and 

?u  /9n.  '^e  hsve  previously  denoted  Kirchhoff 's  solutions  "by  K  and  Z  ,  where  the 
0  as 

su"bscripts  refer  to  diffraction  "by  a  finite  aperture  and  "by  a  finite  screen  res- 
pectively.  It  is  convenient  to  introduce  a  similar  notation  for  the  wave  functions 
derived  from  Rayleigh's  formulas.  "For  ejample,  R  . (P)  will  denote  the  wave 

Six 

function  for  the  diffraction  "by  the  aperture  A  "based  On  Eayleigh's  first  formula 

(l).  For  o"bviou8  reasons,  the  corresponding  modified  Kirchhoff  solutions  will 

"be  termed  "Rayleigh  solutions". 

The  complete  set  of  Eayleigh  solutions  for  the  diffracted  field  in  the  right 

half-space  (z  >  O)  then  "becomes 

.        /    ikr     3u 

(3.3)  B^^(P)  =  «^(P)*i^V    Sf   ^21. 

ikr     9u 


,        /•    ikr     du 

■I  /  a  1^ 

(3.5)         R.,(F)  =  u^(p)  -  i^  u„  Ij  ( V  ^  2: . 


(1.6) 
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TTnlllce  Kirchhoff  8  theory,  the  modified  theory  is  self-consistent^  -' .  The 
reason  for  this  is  that  in  the  modified  theory  it  was  sufficient  to  assume 
■boundary  values  for  either  u  (in  the  case  of  E  )  or  du/Sn  (in  the  case  of  E^). 
In  fact,  all  values  assumed  to  exist  at  z  =  +0,  whether  "behind  the  screen  or 
in  the  aperture,  are  exactly  reproduced  hy  the  Eaylelgh  solutions  when  P  approaches 
the  plane  z  =  0  from  the  right. 

The  analytic  continuation  of  the  Rayleigh  solutions  into  the  illuminated  half- 
space  are  easily  ohtained.  Equations  (3)  and  (5)  remain  valid  for  points  P  to  the 
left  of  the  plane  z  =  0.  On  the  other  hand,  equations  (k)   and  (6)  are  to  he  re- 
placed hy 

,   /•  ikr  du 
(5.7)         H,,(P)  =  u,(P)  -  u„(-P)  -  i/V  ^    ^Z. 


and.  ' 

ikr 
(3.8)    B^^(P)  =  „^(P)  .  u^(.P)  .  i^u„  fc  (V->  ^  ^ 

respectively.  Here  u  (-P)  means  the  valne  assumed  by  u  at  the  reflection  of  P 

0  o 

in  the  plane  z  =  0;  u  (P)  +  u  (-P)  Is  the  zero-order  reflected  field  in  the  sense 

o       0 

of  geometrical  0T)tics. 

LiTce  ICirchhoff 's  solution,  the  Eayleigh  solutions  are  exact  solutions  of 
sRltus  problems .  The  functions  E   and  H  „  Jump  from  2a     on  the  illuminated  face 

S  c,  ct^  0 

of  the  screen  to  zero  on  the  c^arV:  face.  Similarly  the  normal  derivatives  of  E 

8  A. 

and  E  .  Jump  across  the  screen  from  2  9u  /9n  to  zero.  Further,  the  Kirchhoff 

solution  is  Just  the  average  of  the  two  corresponding  Eaylelgh  solutions,  viz., 

K     =  i  (E  ,  +  E  ^)   ,  K  =  i  (E  ,  +  E  ^), 
a   2   al    a2   *  s   2   si    s2  ' 

while  Bahinet's  principle  now  assumes  either  the  form 


\i  +  Kt   =  ^^   (z  >  0),  E   +  E   =  2a„(P)  -  u  (-P)  (z  <  O), 
ai    si    0  al    si     o      o 


or  the  form 


Ko  *  K^   =  ^«  ^2  >  0),  E   +  E   =  2a„(P)  +  u  (-*)  (z  <  0). 
aZ    S/i    0  a2    s2     o       o 
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As  we  mentioned  before,  one  aim  of  the  modified  Kirchhoff  theory  was  to  furnish 
a  method  for  distinguishing  between  the  two  principal  boundary-value  problems  pre- 
viously noted.  Accordingly,  since  3R  -./Bn  =  0  on  the  dark  face  of  the  screen,  R  , 
of  equation  (U)  was  proposed'-  -1  as  an  approximate  expression  for  the  diffracted 
field  behind  the  aperture  in  the  acoustically  rigid  screen  for  the  incident  wave 
u  (i.e.,  R  ^:^  0_),   (^,  and  9*0  sre   defined  in  equations  (9)  and  (10),)  Similarly, 
R  „  was  suggested  for  the  diffracted  field  behind  the  aco-QsticnTly  soft  screen 
(i.e.,  R  p^  ^, ) .  It  is  obvious  that  these  approximations  will  be  accurate  immediate- 
ly behind  the  screen  but  poor  in  the  vicinity  of  the  aperture.  The  approximation 
is  a  complete  failure  if  it  is  extended  to  the  respective   analytic  continuations 
through  the  aperture  into  the  illuminated  space,  because  on  the  lit  face  of  the 

screen  we  have  3R  -,/dn   =  2  du  /9n  and  R  ^  *  2u  ,  In  fact,  the  reflected-field 
al'        0'        a2     o         ' 

terms  in  equations  (7)  and  (8)  suggest  the  opposite  correspondence  between  the 
Rayleigh  solutions  and  the  solutions  of  the  boundary -value  problems.  We  shall 
now  discuss  this  correspondence  in  more  detail. 

Let  u  -,(x,y,z)  denote  the  wave  function  for  the  diffraction  of  the  primary 
field  u  (x,y,z),  incident  from  the  left  (z  <  0),  through  a  finite  aperture  in  a 
perfectly  soft  plane  screen.  Then  '-  -^ 

%(x,y,z)  -  u^(x,y,-z)  +  (jS^{x,y,-z)  (z  <  0) 

{i^lCxjyjz),  (z  >  0) 

where  (2>, ,   defined  for  z  >  0  only,  has   the  following  properties:    (i)  ^^  is  a  solution 
of  the  wave  equation;    (ii)  ^,    =  0  on  the      dark  face  of  the     screen;    (iii)  <^-,    is 

regular  at  infinity  (Sommerfeld's  radiation  condition);    (iv)   b^-.{bz  -  3u  /3z  in 

2 
the  aperture;    (v)  ^     is  uniformly  bounded,   and    |grad  ^-,  |      is  inteq-rable  over  any 

finite  part  of  three-dimensional  space,  including  the  rim  of  the  aperture. 

Let  u  2(x,y,z)   denote   the  corresponding  wave  function  for  an  aperture  in  a 

perfectly  rigid  screen.     Then 


(3.9)         u^ 
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\(^,7*z)   +  u^(x,y,-z)  -  ^^{x^y,-z)  (z  <  o) 


(-^.10)         Tl.^  = 


a? 

'  0,(x,y,z),  (z  >  0) 


where  f)     is  also  defined  for  t;  >  0  and  has  similar  properties  as  ^.   except  that 
(ii)  and  (iv)  should  "be  replaced  "by:  (ii')  ^2^^^   =  0  on  the  screen,  and  (iv») 
02  =  11  in  the  aperture.  Existence  theorems  and  questions  of  uniqueness  will 
not  be  discussed  for  the  time  "being.  Let  it  suffice  to  mention  that  the  pro- 
perty (v)  ensures  that  no  energy  is  radiated  hy  the  singularities  at  the  rim. 
It  is  difficult,  if  not  impossible,  to  determine  the  functions  0-  and  0_ 
for  an  aperture  of  arbitrary  shape.  The  trouble  is  that  either  0^  or  0  solves 
a  mixed  boun dp ry- value  problem:  0  and  90 /Sn  are  given  on  mutually  complementpry 
parts  of  the  plane  z  =  0  (see  (ii)  and  (iv)).  However,  by  virtue  of  (ii)  and  (ii')' 
and  Rayleigh's  formula,  we  hnve  for  any  aperture  A  the  following  relations: 
-   /.    ;.   ^ikr  -   /»50o  -i^ 

!^ov'  if  we  assume  thpt  the  unknown  vs. lues  of  0.  and  90-/Bn  in  the  aoerture  Bay  be 
replaced  by  the  respective  unperturbed  values  of  the  incident  wRve,  we  find  that 
^T  ~  ^-ot  ''-ft*  R  .  (z  >  O).   If  this  Is  substituted  into  equations  (9)  and  (10 ) 

J.      8.C     ^      SIX 

we  ©"btftln 

^al^\2^^>°^'    ^al*^o-\2'   ^^  <  °^ 
(3.12) 

\2-\l^^>^^'         ^a2«2u^-\l'   (^'^°^- 
For  z  >  0  the  approximation  (12)  is  identical  with  that  discussed  previously. 
It  should  be  noted  that  the  approximate  solritions  do  satisfy  the  correct  boundary 
conditions  at  the  screen.  However,  they  are  not  analytic  functions:  either  the 
normal  derivative  or  the  approximate  solution  itself  is  discontinuous  in  the 
aperture. 

In  deriving  the  approximations  (l2)  the  properties  (ii)  and  (ii*)  have  been 
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used.  Tet  almost'  eq>»lly  characteristic  for  planar  boundery- value  problems  are 
the  properties  (iv)  and  (iv')  .  The  latter  eacpress  theit  for  problem  I  (u  =  0  on 
the  screen)  9u/?^n  is  unperturbed  in  the  aperture.  Starting  from  this  point  of 
view  we  can  derive  a  different  set  of  approximations'-  -■  .  From  Eayleigh's 
formula  (l)  and  property  (iv)  it  follows  that 

,  .    ikr  3u         ,    f  J-"^     ^01    __ 

Hnfortunptely,  the  values  of  c^()-/9n  are  not  known  on  the  infinite  screen  S,   If 
we  assume  that  they  are  approximately  equal  to  zero,  we  arrive  at  0.  s?  ^oi'  -^^ 
a  slmilnr  reasoning  in  the  case  of  problem  II,  we  ^et  0^  *  -^  p«   -^^  ^^  J^o*  diffi- 
cult to  verify  that  this  ultimately  gives 

t 

(MM        ^al«  \l  •   ^a2=^\2 

everywhere  in  free  space.  The  approximate  solutions  are  now  analytic  end  they 

produce  the  correct  (imperturbed)  values  of  du  n/Sn  and  u  g  in^  the  aperture,  but 

they  violate  the  correct  boundpry  conditions  at  the  screen.   Insofar  as  an  accurate 

approximetion  to  the  field  is  more  important  in  the  aperture  than  in  the  vicinity 

of  the  screen,  the  approximations  (l^)  seem  preferable  to  the  opposite  relations 

(12). 

An  alternative  way  of  sho\irin#'  the  close  relption  between  the  original  and 
modified  Kirchhoff  solutions  is   the  following.      Consider,   for  example,   the 
diffracted  field  behind  an  aperture.     Noting  that  a/Bn  =  S/Sz '  =  -9/3z  for  any 
function  of  r,  we  hpve  from  equation    (^)  that 

(''.lO  2E^(P)  =  f(P)  +  ^    g(P), 

where  f(P)  and  g(P)  are  both  even  functions  of  z,  viz., 

-    t    ikr  ?u  ,   /  ikr 

(..16)        f  (P)  .  -  ^  _^  V  ST  *  E  .  .(F)  =  -  hi'—    V  JI- 

Comparison  with  (^)  and  (6)  shows  that 
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(3.17)        \l(P)  =  f(P),   E^2^^^  "  h    ^^^^- 

Consequently,  the  functions  R  -.  (P)  8ind  ^2^^^  ®^®  simply  twice  the  even  and  odd 
parts  respectively  of  K  (P). 

In  addition,  the  function  f (P)  is  equal  to  the  velocity  potential  of  a  mem- 
brane vibrating  in  a  rigid  haffle  with  velocity  distribution  -9u  /hn;   the  same 
holds  for  g(P)  and  -u  .  Methods  and  n\3jnerical  results  of  the  theory  of  acoustic 

0 

radiation  are  therefore  valuable  in  diffraction  theory  also.  Vsrious  suthors 

ri2i 

have  used  the  modified  Kirchhoff  theory  in  one  way  or  another.  Bremmer"-  -*  applied 
Reyleigh's  second  formala  to  the  diffraction  theory  of  Gfl.us8ian  optics.  Various 
mBthemstinal  aspects  of  equation  (2)  were  discussed  by  Luneberg'-  ■-' ,  and  Scheffers'-  -^ 
emphasized  the  usefulness  of  this  6quation(in  the  Fourier  form)forthe  calculation  of 
Fraunhofer  rattems.  Durand^  -'  applied  the  Fooorier  equivalent  of  (6)  to  a  circular 
aperture  end  to  p  helf-plane  for  the  case  of  a  plane  wsve  at  normal  incidence. 
Spence-  -  compared  R  ^  for  a  circular  aperttire  (plane  wave  at  normal  incidence) 

Si  J. 

with  the  corresponding  exact  solution  of  the  boiindary- value  problem  II.  Experi- 
mental results  on  the  diffraction  of  sound  around  a  circulardisk  were  discussed 
in  connection  with  the  Kirchhoff  approximations  by  Primakoff ,  Klein,  Keller  and 
Cars  tens  en'-  -'  , 

In  concluding  this  section,  we  shall  briefly  discuss  the  Kirchhoff  solutions 

for  the  diffraction  of  a  plane  wave  normally  incident  on  a  circular  aperture.   Let 

ikz 
a  denote  the  radius  of  the  aperture,  and  let  the  incident  wave  u  =  e    impinge 

from  the  left.  Choosing  the  origin  of  coordinates  at  the  center  of  the  aperture, 

we  have  in  the  shadow  region 

(..18)  ^   ^  /e^  ,V-  R   -   ^ 

where  TT  is  the  velocity  potential  of  Rayleigh's  piston  for  unit  velocity  dis- 
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tritcutlon. 

The  integral   is   p-asy  to  evaltipte   if  P  lies  on  the  z  axis.     The  result   is 

f^  .,.\  -D  i'^z         ikyz^+a^ 

(^.19)  \-i  =  ®         -  e       '  , 

al 


'z 


-4.a* 


/oo■.^  ir  i^       l^ikYz^f?       1         «  ^ik /z^+a 

(^.21)  '^^  =  e         -  9  *  -  ?      .  e 

*  ^  /  2     2 

yz  +8^ 

It  should  bf!  noticed  that  these  ea^jressions  are  equally  vAlid  if  P  is  on  the 
negative  z  axis.   If  a  tends  to  Infinity,  R  ^  and  K  do  not  reduce  to  the  inci- 
dent field  e  ^,  unless  the  medium  is  assujned  to  he  slightly  absorbing  (i.e., 
imaginary  part  of  k  is  positive). 

The  respective  PraTinhofer  patterns  are  also  easy  to  calculate.  Let  p,  6 
denote  spherical  coordinates  with  the  positive  z  axis  as  polar  axis.  Then  at 
ler^e  distances  from  the  aperture 

U-^A(e)  p-^  e^'^P  , 

where  A(o)  is  the  amplitude  of  the  spherical  wave.  We  find 

J^ (kasine) 
(..22)         (i/a)A^^  =  %ine    ' 

J  (kasinQ) 


» 


J- (kasinO) 
C^.PM         (i/a)A^=   2tan(Q/2)   » 


where  we  use  an  obvious  notption  for  the  amplitude,  and  where  0  <  6  <  —  ,  and 
J^  is  a  Bessel  function. 

The  amount  of  energy  transmitted  through  the  aioerture  can  be  computed  by 
integrating  |a|   over  half  of  the  unit  sphere.   It  is  convenient  to  introduce 


the  transmlflBJon  coefficient,  which  1b  the  ratio  of  the  trsnamltted  energy  to 
the  energy  Incident  on  the  apertuT-e  in  the  sense  of  geometrical  ontics.   In 
the  prohlem  nnder  discussion  this  coefficient  is 

Tr/2 
(3.25)  T  =  -^   /lAJ^dA  =  -%  /     |A(«)|28in«de  . 

The  relevant  expressions  for  the  three  different  cases  mentioned  ahove  are 

J,  (2Tca) 


(?.26)  T,  =  1  -  ^ 


1        ka    • 


J.  ( 2kB. )    ,    2ka 
(-^.27)  T^  =  l^.^__.^  A  j^(Odt. 

''o 

(?.28)  ■   T J  =  1  -  I  C  jf(ka)  +  J^O<a)  +  ^  /    J^(t)dt] 


Equation  (26)  Is  a  classical  result  ohtained  by  lord  Rayleig'h. 
""'or  very  small  values  of  'l<a.  we  heve 

^.29)        T^-  |(ka)2  ,  T^-  |(i^)2  .  ^K-ir  ('^)^- 

These  -^.lues  of  T  are  in  complete  disegreement  with  the  results  for  the  exact 
>>oundary-value  prohlems.  This  is  not  surprising,  since  the  Kirchhoff  approxima- 
tion holds  for  small  wavelengths.  For  very  large  values  of  ka  we  hs.ve 


m    -    1_     sin(2ka)-TT/It) 
^^^  "  r-  ,.     v3/2   • 


Vn       (kB.)- 


(3.30)  ^2-^-i^-^^  cos(2ka-n/l;) 

^       ^        Zfi  (ka)5/2 


I 


m        ,   fi  +  ll  i_   _1_  3ln(2ka)-TT/li) 
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IV,  BraimT)e>:'g  Modification  of  the  Klrchhoff  Theory. 

An  attempt  to  Improve  the  modified  Klrchhoff  8ol-:itlon  was  made  "by 
Bratmhek''  ■■  who  ohserved  that  the  Klrchhoff  solution  does  not  constitute 
the  main  term  (in  a  aeries  of  powers  of  l/ka)  of  the  exact  solution.  How- 
ever, "before  entering  into  a  discussion  of  Braun'bek's  theory,  we  shall 
first  discuss  the  solution  constructed  "by  Sommerfeld  for  the  diffraction 
of  a  plane  wave  by  a  half-plane: 


£=  Ae 


-ncpcos(e-e^)   Vacpcosfe-ej/zl  ^^2 


ya^< 


dx 


-ikpcosCdfO    )    yapco8[(©+©  y23  2 

+  Be  °    /  e  dr    , 

on  which  Braunbek's  theory  is  "based, 

A  simple  derivation  of  Sommerfeld 's  f omuls s  is  implicit  in  an  inter- 

ri9i 

esting  paper  by  J.  Brillouin'-   ;  this  proof  we  now  present. 

Consider  a  plane  wave  incident  on  a  screen  which  cuts  the  plane  of 
drawime  in  the  Mpver  half  of  the  y  axis  (fig.  1,  direction  of  incidence 
normal  to  edge  of  •screen). 


Screen 


ng.  1 

(The  restriction  to  normal  incidence  causes  no  lack  of  generality  since 
the  proof  can  be  generalized  to  arbitrary  oblique  incidence.)   Introduce 
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new  "Bemi-para'bolic"  coordinates  u,  v  such  that 


(J^.l) 


▼  =  7 


n  =  -7  +  Y7^+B^ 


7  »  T 

z 


^S/^' 


+Zar 


(i.e.,  one  of  the  psraholic  coordinates  is  rejected  in  favor  of  a  rectangnlar 
coordinate),  which  can  he  written  in  terms  of  polar  coordinates  as  well: 

▼  =  -  p  cos  9 


(^.2) 


2  © 
u  =  2p  cos  2  , 


where  0  <  p  «<  oo  and  0  <  «  ^  rr. 
Now,  from  (l)  we  see  that 


(^4.3) 


a   z  9 

^■■1  SE   ■!■  ■■     I  »■     • 

az  u+v  a-a  * 

,2 


37   av 


n  a_ 
u+v  au  * 


With  the  ttse  of  (3)t   A  ==   9  +  ^"^  hecomes 

ay^  az"^ 


(^.i*) 


ar^   ^ 


au 


s_  _  a' 
2  ~  avau 


>  + 


v+u  au  * 


The  wave  equation   ^J+k  5=0  is  not  separable  in  these  coordinates; 
however,  we  follow  the  same  procedure  as  we  would  if  it  were  sepairahle  and 
assume  that  £  =  F(v)G(u).  The  wave  eqti&tion  then  hecomes 


(i».5)     Glj"+k^  +  ^ 


\r*-v 


P'  &"+  ^  G«f  - 


r«G« 


=  0 


which  can  he  satisfied  if  we   take  T  and  G  such  that  the  relations 

(i'.6a)  ?•»  +  k^  =  0 

(ij.6h)  y^G«+:^G«^    >P«G'  =  0 


?/  G«+  i-  G«l     -  P'G'  = 


are  satisfied. 

T'rom  <6a)   it  follows  that 

ikv 


(^.7) 


y  «  a  e' 
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and  when  this  result  is  stihstltuted  Into  (6"b),  the  latter  'becomes 

It  is  seen  from  (8)  that 

G«  =  ^ 


and  finally  that 

u   ikt 

(2J.9)  G  =  P/  7---  it  . 

•^  Ut 


Using  the  results  of  (?)  and  (9)  the  solution  5=  F(t)G(u)  becomes 


(^.10)        FG  =  (const.)  e^ J 


u  ikt 


\A 


Transforming  to  polar  coordinates,  we  obtain 

2pcos^(e/2) 

(a.ll)        TO  =  (const.)  e-^'^P^^^V        S_^  dt  , 


>^ 


and  moreover,   If  we  let  kt  =    Z   ,    (ll)  "becomes 


41    -  -fl    y2^cos(©/2)     ,  _2 
{U,12)  rG=   (const.  )e-^'^P^°''Y  e^  ^     dT   . 

If  we  had  started  with  Q  +  ^     instead  of  Q   (this  just  means  a  rotation 

~  0 

of  the  coordinate  system)  we  weald  have  obtained  two  other  solutions,  and 
by  taking  a  linear  combination  of  these  two,  we  would  have  as  our  final 

-lkpcos(G-9  )  jf^cos-^ 
(it. 13)        FG  «  Ae  °  y  e^  ^   dT 

o 

-ikpcos«He)  y^cosi^ 
^■Be         °  y  e^         dC- 

o 
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which  Is  exactly  the  form  of  Somraerfeld'a  solution  and  hence  verifies  it, 
(The  "boimdary  conditions  can  he  met  by  an  adjustment  of  the  constants  A 
and  B.) 

We  shall  now  treat  the  two  special  cases  referred  to  in  Section  III 
(i.e.,  ft.   =  0   on  the  screen,  and  hft^^/hn  =  0  on  the  screen).  However,  we 
shall  use  henceforth  a  slightly  different  system  of  coordinates;  see  fig, 
2. 

Incident     Wave 


->Z 


Fig.  2 


Defining  ^-  and  f)^   to  he 


^l] 


=  e 


ikpcos(e-«  )  ^-i"/i4  8^  ^^2    _  ikpcos(d+e^)^-i7ii  8   ^^2 

•  /   e    dT  +  e  /   e    dT 


v;;  < 


v^  ->. 


Q-Q  *4-©  lkpcos(e-9  ) 

where  s  =  V2ko  sin(  g  °),  S-  =  -/2kp  sin(  ^  °),  and  where  u  =  e         ° 

is  the  incident  wave,  it  can  he  verified  directly  that  >5^  and  30-/9n  hoth 
vanish  on  the  screen. 


There  are  other  fijnctions  which  satisfy  the  wave  equation  and,  at  the 

same  time,  the  boundary  conditions.  Thus,  for  example,  the  functions 
9  jfi,/Sz  ,  3  P-tl^^    t  etc.,  which  obviously  satisfy  the  wave  equation,  also 
vanish  on  the  screen'-   .  However,  these  functions  are  too  singular  at 
the  edge  and  hence  are  not  admissihle  solutions;  fi^  and  0p  of  equation 
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(lli)  are  admlesible'-^^J  . 

For  the  purpose  of  discussing  Sommerf eld's  theory,   two  electromagnetic 
fields  can  "be  constructed  from  the  solutions  of  the  scalar  waye  equation 
previously  discussed.     The  first  field  is 

(^.15a)  ^11    =   (1*1,0,0), 

which  vanishes  on  the  screen,  and  for  which  H,  determined  "by  Maxwell's 
equations,  is  given  hy 

and  the  second  field  is 


(li.l-?^)        H|  j  =  (162.0,0), 

the  normal  derivative  of  which  vanishes  on  the  screen,  and  for  which 

30,  -30J, 

The  explicit  expressions  for  all  field  components  for  the  Sommerfeld 
half -plane  prohlem  are 

\=^1 

Ti      ^    •  Q        1  fl        1  O  r^  i(kp+TTM) 

H  =  ?5,co86  -cos  ig  Sin  ^9  V^  ^H^p^/^) 
z   ^1    o     2  0     2   'Trkp 

in  the  case  of  incident  field  polarised  parallel  to  the  edge,  and 

for  the  incident  magnetic  vector  parallel  to  the  edge.  For  the  defini- 
tion of  coordinates,  see  Pig.  2,  not  ^ig,  1. 
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Braunbelc'8  theory  makes  use  of  the  values  of  the  scalar  fields  and  their  normal 
derive  tires  on  the  x.y  plane  (i.e.,  the  plane  determined  hy  the  screen). 
The  non-vanishing  derivatives  of  i^   and  ^^  necessary  to  compute  the  field 


quantities  are 


3z 


(iksine^)f52  4.yf  e-^  ^""^o^lej 


cos*  d  - 


ikp 


rp 


dz 

(^.16) 

3|«1 


(iksinO  )«5, 
0   1 


^— e  '   8in-e^^8in|©  ^ 


Ikp 


V? 


g— =  (ikcose^)^^  - 


llT  ®  '  cos  -  e^  ( 


1   e^'^ 
sin^  Q   2 

'      Vp    J 


sin  -  9  \co8-  Q 


To  calculate  the  fields  on  the  x,y  plane,  we  let  e  =  0  and  e  =  tt. 


We 


confine  ourselves   to   the  case  of  normal   incidence  and  hence   tsfce  ©     =  - 

0       2 
This  ^ives  u8 


(^■.17) 


and 


(^.18) 


©  = 


^1 


=     0 


=     0 


'< 


^2  ~V^ 


_2     -i^/W-^^        iT^ 


az 


=  ik 


©    =    TT 


/ 


h    =1 


VTT  ^k^ 


TT/i  00  .         2 

e         dr 


?^  =  ii.  -  m  .-!■■/'• 


V    9z 


>fn 


/        e  dx     K  ,^ 


^ 


/T 
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We  shall  now  discuss  Bra-unhek's  attempt  to  improve  the  modified 
ICirchhoff  solution  for  a  circular  aperture.  As  was  mentioned  in  the  "be- 
ginning of  this  section  Breunhek  observed  that  Kirchhoff 's  solution  does 
not  constitute  the  main  term  of  the  exact  solution.  To  obtain  the  correct 

main  term,  we  must  estim^-te  the  effect  of  the  second  integral  of  (3.13) 

901 
on  0  .   (The  problem  for  ^g  ^^  analogous.)  Braunbek  replaced  g^  on  the 

screen  by  the  value  derived  from  Sommerf eld's  theory  for  the  half-plane, 
as  if  S  hpd  p  locslly  straight  edge.  This  is  a  plausible  assumption,  be- 
cause hi^/bn   la  ejcpected  to  be  rapidly  decreasing  from  the  edge  over  a 

distance  of  a  few  wavelengths^conf irmed  experimentally  in  the  case  of  f)^ 
by  Severln  and  Starke^  y^   and  the  wavelengths  here  considered  are  small 
with  respect  to  the  radius  of  the  aperture.  For  the  circular  aperture 
and  a  normally  incident  plane  wave,  Braunbek's  approximation  becomes 

(..19)  i,«\  =  -  i^/^^E-  i^^  [l(k!)-.fel)|  dE  . 

where    t    is   the  distance  to  the  rim  and 

/.        N     T/    N        2      -i"M    r°°   iT^  ,    N       e^  /^     ix 

(i^.20)    i(x)  =  -!:  e  ^  '     /_  e       dr  ;  mW  =  -——  e      . 

The  evaluation  of  B  on  the  axis  of  the  aperture  is  comparatively  simple. 
Integration  by  parts  and  some  trivial  transformations  yield 

(ii.2l)       B^(0,0,z) 


-  e*"*^^ 


•r-4-va-*-kyz  +a 


®        d-c. 


J+2k\/z^+a^'j  Yz 

The  integral  is  elementary  if  z  =  0.     At  the  center  of  the  aperture  Breunbek's 
function  is  exactly  equal  to 

(i».22)  B^(o,o.o)  =  1  -V2  e^^. 

When  z  Is  greater  than  zero,   the  integral  in   (21 )  can  be  expaartLed  -oymptotically. 
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In  a  first  approximation,  the  square  root  in  the  integrand  may  be  replaced 
■by  its  value  of  the  integrand  at  the  lower  limit  of  integration.  Then  ex- 
cept for  terms  of  order  l/ks.. 


^V^ 


2._2. 


(il.23)        B^(0.0.O=  e^^^  -  e^^V^  "^  C 1  +  -J^l^^^   . 

\z   +a 

which  is  Braunbek's  result    (obtained  somewhat  differently).     His  computations 

showed  that   (22)  and   (23)  are  in  excellent  agreement  with  Meixner  and  Tritze's  "-^^-^ 

exact  values.     The  corresponding  diffraction  pattern  was  also  evaluated  by 

Braunbek. 

In  studying  Brannbek's  paper,  the  author  encounters  some  difficulties 

in  connection  with  a  second  type  of  approximation  to  ^. ,   In  starting  from 

the  equation 

given  in  Section  III,  Braunbek  replaces  ii.    In  the  aperture  by  the  value 
1  -  5(k|)  derived  from  Sommerfeld's  theory.  Then 

Bratinbek  claims  that,  except  for  terms  of  order  l/ka,  the  functions  B  and 
B*  are  identical  on  the  axis  of  the  aperture.  The  present  author  believes 
this  to  be  incorrect,  because  from  the  values  at  the  center  of  the  aperture 
we  can  see,  without  any  calculation  at  all,  that 

B*  (0,0,0)  =  1  -  5(ka)  -  1  -  7—  e^^  , 
which  differs  from  (22). 
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T,     Tarlatlonel  Formulation  of  Scalar  Diffraction  ProblemB. 

A  vprlatlonal  formulation  of  planar  diffraction  pro'blems,  which  permits 
accurate  numerical  evaluation  of  the  diffracted  amplitude  and  the  trans- 
mlasion  cro«3-section  for  a  wide  ren^e  of  frequencies,  was  given  "by  Levlne 

r??i 

and  Schwln^er'-  ■■.  They  illustrated  the  utility  of  the  variational  method 

y>j  applying  It  to  the  circular  aperture  for  a  normally  incident  plane  wave. 

r2i4l 
The  analysis  wss  criticised  by  Copson"-  -'  since  many  of  the  integrals  in- 
volved appeared  to  diverge.   Copson,  however,  in  deriving  what  he  calls 
"Levlne  and  Schwinger's  variational  principle"  in  a  mathematically  sound 
way,  confined  himself  to  the  prohlem  for  0  ,  while  his  criticism  concerns 

ths  t  for  0. ,  ■^''ortunatfcly  the  divergent  integrals  that  occur  in  Levine 
and  Schwinger's  paper  are  easy  to  eliminate  without  affecting  the  numerical 
results.  However  before  we  treat  these  pro'blems  let  us  first  discuss  the 
integral  formulation  of  scalar  diffraction  pro'blems  with  application  to 
small  apertures. 

As  is  seen  from  the  second  eqiiation  of  (3.1l)»  the  wave  function  ^_ 
is  uniquely  determined  in  space  'by  the  values  of  its  nonnal  derivative  in 
the  aperture.   Let  the  unVrnown  values  of  S0  /dn  in  the  aperture  'be  denoted 

hy  f(x,y).  Recalling  that  )d„  =  u  in  the  aperture,  we  find  the  integral 

[25] 
equation'-  ■' 

(5.1a)       /  f(x',y')G(s:.x',y,y')dx'dy'  =  -u^(x.y,o), 
^  0 

where   the  kernel  G  is   sjrmTiietrlc  and  singular,  viz. 

G  =    (V2TT8)e*^  ,        8^  =    (x-x')^  +    (y-y')^   . 

and  where  (x,y,o)  Is  any  point  of  the  aperture. 

Similarly,  the  first  equation  of  (3.11)  shows  that  0^  is  uniquely 

determined  'by  its  value  in  the  aperture: 


1   dz 


-  2^^  - 

Tn  thlfl  cane  we  h«v©  the  further  condition  thet  3)6. /9z  =  Su.  /5z  in  the 
aperture.  Since  in  eonption  (l)  the  term  in  "brackets  is  a  solution  \y, 
fi^-Jt   of  the  wave  equation,  we  have 

oz  ox         oy 


80   that  the  following  relation  is  o'btained 


M. 


9u  (x.y.o)     ^   ^2    2,2  _,  , 
(5.2)     —5-^;; =  &""  "^  ~  +  ^-rH  /  )«,  (x '.y  «,o  )G(x.i',y,y '  )dx'dy', 

where  (x,y,o)  is  in  the  aperture.  This  relation  is  not  a  pure  integral 
eq_u»tion;  it  is  a  differential-integral  equation.   It  should  he  noted  that 

the  differential  operator  may  not  be  applied  to  6  under  the  integration  sign 

[271 
since  the  resulting  kernel  would  not  he  integrahle.  Maue^  -*  gave  an  equi- 
valent form  of  equation  (2),  namely, 

(5.?)   ^    =  ^^    /  ^i^^  ^  "/  ^^^*^  ^1   '  ^^^     '^'^^  ^  • 
where  "both  gradients:  are  to  he  taken  with  respect  to  the  coordinates  of 
integration  x',  y'.  Equation  (3)  follows  from  (2)  by  a  process  of  differ- 
entiation and  integration  hy  parts,  and  use  of  the  condition  0-  =  0  at 
the  edge  of  the  aperture.  The  second  integral  in  (3)  is  a  Cauchy's  - 
principle  value  (small  circle  around  x,y,o  of  radius  €  ->  O).  A  second 
integration  hy  parts  is  impossible  because  of  the  singularities  at  the 

edee. 

Only  in  a  few  simple  cases  can  the  differential-integral  equation  (2) 

be  transformed  into  a  pure  integral  equation.   If  the  incident  field  is  a 

plane  wave,  u  =  exp[j.k(ax+Py+Vz2  ,  "V  =  (l-a  -P  )  '  ,  we  have 

!l4:^  .  ^,2  ,  iL,  iL^  (,/,,)  3ik(a^Py)  ^ 
^*  ax^   3y2 

so  that  (2)  becomes 
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(5.M   /0^(x'.y'.o)G(x.x'.y.y')  =  ^  ^^^(c^^y)  ^   ^(^.^y)  , 

where  /    is  a  solution  of  the  twb-dimensionsl  wave  equation  in  the  aperture. 
The  ftinction  Y  can  be  uniquely  determined  (except  for  a  constant  multiplier) 
for  a  normally  incident  vrave  (a  =  P  =  o)  in  the  eases  of  the  circular  aper- 
tiire'-  •'•  and  the  infinite  slit'"  ■',  The  resulting  equation  can  he  trans- 
formed into  a  pure  integral  equation  with  a  non-symmetric  kernel. 

Let  us  now  return  to  the  discussion  of  the  variational  formulation  of 

planar  diffraction  prohlems  as  given  hy  Levine  and  Schwinger.  As  was  mention- 
ed "before,  the  divergent  integrals  that  occur  in  their  paper  are  easy  to 
eliminate  without  affecting  the  numerical  results.   In  what  follows,  the 
notation  is  stiggested  hy  that  of  Levine  and  Schwinger. 

Let  A  (n,n')  denote  the  amplitude  of  the  diffracted  wave,  where  n  and 
n'  are  unit  vectors  in  the  direction  of  observation  and  of  propagation  of 
the  Incident  plane  wave  respectively.  Further,  let  9   and  ©'  "be  the  angles 
"between  the  positive  z  axis  and  these  unit  vectors.  From  equation  (l)  it 
^'ollow8  that 

(5.5)  A^(n.n')  =  -  (ik/2TT)co3  6^  0^,  (pOe"^^"^ 'dS  • , 

where  6   t(p')  is  the  value  of  f>,   in  the  aperture.  The  integral  equation 
n  1 

in  Levine  and  Schwinger*8  paper  [their  equation  (2.911  contains  a  non- 
integrpble  kernel  and  should  be  interpreted  in  the  sense  of  one  of  the 
differential- integral  eq-uations  (2)  or  (3).  Let  us  choose  Maue's  equation 
(1),  Then 

(5.6)  2TT  ik  cose'e^^'P  =  k^  /0^,(p')G(p,p')dS'  -  /r'^^^,  (p ')  7'G(p.p')dS  • . 


where  p  is  in  the  aperture  and 
(5.7)         (iip.p')  = 


ei^^lp-p'l 


Ip-P'l 
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If  we  nittltlply  throneh  in  eqiiation  (6)  by  0  „(p)  (where  0  ,t(r)  is  the 
solution  for  a  plane  wave  in  the  direction  n"),  and  integrate  over  the 
aperture,  there  results 

(5.8)  2TTiVcose'  /^«)^„(p)e^^'PdS  =  k^  0^„  (p)G(p,p«)0j^,  (p*  )dSdS  ' 

-  /^^„(p)^')!5^,(p')G(p,p')dSdS». 

The  last  integral  appears  after  an  integration  hy  parts;  the  integrated 
term  drops  out  because  0  „,  like  jd  ,,  vanishes  at  the  edge  of  the  aperture. 
The  right-hand  side  of  eqijation  (8)  is  symmetrical  in  n'  and  n"  and  con- 
sequently so  is  the  left-hand  side.   If  we  divide  (8)  by  this  left-hand 
side  and  by  a  similar  term  in  which  n'  and  n"  ere  interchanged,  and  then 
use  (5)  and  Invert,  we  obtain 

(5.9)  A^(n",n')  =  A^(-n«.-n") 

cos«'cos9"/0^,(p)e-^^"PdSy^_^„(p)e^^'PdS 

'  /[^n'^P^^-n"^P'^-^'^''^n'^P^'"^-n"^'^'^j^^P'P'^'^^<^^'  ^ 

It  can  be  proved  thpt  the  equation  (9)  is  stationary  with  respect  to  small 
Independent  variations  of  i>    ,  and  0  „   about  their  correct  values;  those 
variations  which  do  not  violate  the  condition  0.  =  0  at  the  edge  are 
admissible. 

The  stationary  character  of  the  expression  (9)  is  of  importance  for 
approTi'Tate  calculations.   In  fact,  a  judicious  choice  of  aperture  distri- 
bute-"ns  in  equation  (9)  may  result  in  a  reasonably  correct  value  for  A 
without  the  necessity  of  solving  the  original  differential-integral  equation. 
Scale  factors  are  of  no  acccont  since  (9)  is  homogeneous  of  degree  zero  in 
0,  The  same  remarks  apply  to  the  plane-wave  transmission  cross-section  <3^ 
of  any  aperture  in  a  plane  screen  (perfectly  rigid  or  soft)  whiebj accord- 
ing to  Levine  and  Schwinger ^-^■' ,  is  related  to  the  amplitude  A  of  the 

spherical  wave  at  large  distances  behind  the  aperture  in  the  direction  of 
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th©  Incident  wave  "by 


2TT 


^     =  .  ^     In,  A. 


"P'rom  equ-'>tlon   (10 )  It  follows  that 


(5.11)   0\(n)  =  -  |^co8^e  Im 


y0^(p)e-'^PdSy0_^(p)e^^PdS 

y|0^(p)j«_^(p')-k"^V()^(p)V'0_^(p)]G(p.p>)dSdS' 


Levine  and  Schwinger  discussed  the  limiting  form  of  this  eq.\iation  for 
low  and  high  frequencies.   In  the  static  case  as  k  approaches  zero 


(5.12)  G  = 


ll^lp-p'l 


IP-PI 


-H-+ik-|-|p-p'l-^lp-p't'-o(k^ 
Up-p'1     ^ 


Usin^  equation  (12),  we  obtain  the  following  for  the  denominator  of  (ll): 
(5.13) 


I P-P ' I  1 P-P  '  I 


+  ik^(0  0  +  z!p-p'I^V0.v'0  )   -KjCk^) 


n  ^-n 


dSdS • . 


n^-n  6l^-P 
At  this  point  we  shall  prove  the  theorem 

(5.1^)  /&x6.7jdJi'i7*S^i^U,y)'\/'()^(T\y')7ir,T',7,7')  = 

which  is  necftssnry  for  the  fiirther  calctilation  of  expression    (I3).     We 
assnme   that  0^  and  0_  are  ar'bitrary  functions  defined  in  the  region  S  and 

equpl  to  zero  on  the  hotindary  of  S.      The  left-hand  side  of  equation  (l3) 

can  he  written  in  the  form 

(5.15)         Mxdy  V^^U,7)   •  /dx'dy'  <  V  {  02^}   -  02  V^  [   • 

However,  it  follows  from  partial  integration  and  the  condition  that  0=0 
on  the  boundary  of  S  that 
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/dx'd7'  V  {^/\  V        02^(direction  cosine)  =  0. 


Bouadp  ry 
Hence,   expression   (15)  'becomea 


-  /  dxdy     ^0^(x,y)   •   /dx'dy'02(3r'.y*)     ^'^t 


which  can  "he  written  as 

(5.16)  -  /dx'dy'jJgC^c'.y')   /drfy  ^  0^(x,y)     V'T. 

TTsln^  a  simllsr  procednre  on  the   integral  of     ^0^  (x,y)     T7'F,   the  left  side 
of  equfltlon   (l'*)  can  he  shown  to  he  eqnal  to 

+  /dx'dy'   i^{T\7^)J  diTidj  03^(x,y)     ^   ^F. 

If  we  assume  that  F  is  a  function  of  an  argument  of  the  form  (x-x')  +(y-y') 

then   ^7T7»  =  _y  =  _  7«,  and  thus  the  theorem  is  proved.  Therefore  if  F 
is  of  the  form  |p-p'|   then  we  find,  making  use  of  the  theorem  (l**),  that 

Now,  to  continue  our  discussion  of  the  limiting  form  of  eauation  (11 ), 
we  can  see  that  the  expression  for  the  denominator  becomes 

(5.17)   [:-/!?i^*  ^2/^:1^^  i  ik^ySj  ^  *  o(k^: 

/       |p>p.l  J    |p.p.|      3  y   n-n 

The  numerator  of  equation   (ll)  hecomes 

(5.l8)[yi5^?5_^  +  lien  ■/(o-p')iJ_^-^^^/t>J_^  &  .   (p-p'D^JdSctS' 

Howerer,  for  small  apertures  (with  respect  to  the  wavelength)  0  =  {5   =  {5 
since  in  this  cnse  0  is  not  dependent  on  a  special  n.  Therefore  the  lepding 
term  of  the  transmission  cross-section  is  given  "^y 
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2n,k       2„   n/^(p)dCl 
3 


(5.19)   (3'^(n)-'Y  k  cos'e  =^-^-^- as  k 


^         I 1 1 


IP-P 
which  is  an  illustration  of  Eaylei^h's  \       law  for  email  scatterers, 

r3ii  r32l 

The  analogous  problem  for  (jf^t   discussed  by  Copson'-  •"  and  by  Levine  ■-  ■■ , 
is  easier  since  it  is  based  on  the  pure  integral  equation  (la).  From  this 
equation  it  follows  that 

(5.20)  yY^,(p<)G(p,p')dS<  =  -2iie^^'P, 

where  Y  ,(p')  is  the  value  of  3Q)„/3z  in  the  aperture  due  to  a  plane  incident 
n'  ^  - 

wave  traveling  in  the  direction  n' .  The  amplitude  can  be  found  from  the  second 
equation  (3.11)! 

(5.21)  A2(n,n')  -  -  i^/l^,  (p- )e-^^P'dS', 
and  the  analogue  of  (9)  becomes 

>  ,(p)e-^^""PdS  A  „(p)e^'^'PdS 

(5.22)  iU(n",n')  =  A,(-nS-n")  =  -^ — , 

/Tn,(p)G(p,p')Y_^„(p')dSdS< 

which  expression,  again,  is  stationary  with  respect  to  small  variations  of  Y  , 

and  Y  _  about  their  correct  values.  In  this  case  the  variations  are  not  res- 

-n* 

tricted  by  a  condition  at  the  edgej   the  correct  aperture  fields  Y  are  infinite 
there. 

The  corresponding  formula  for  the  transmission  cross-section  is 

2.  A(P)e-^'"^<is/l     (p)ei^PdS 

(5.23)  a',(n)   =  -  ^  Im  :^^^ :^l-=2 , 

yY^(p)G(p,p')Y_^(p')d8dS' 

while  the  leading  term  in  the  static  case  is  given  by 
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[>(p)ds]^ 

(5.2U)  &^{n)^2ji    (k-».e). 

\y  ip-p'i       J 

2 

This  leading  tenn,  it  may  be  noted,  is  equal  to  2tiC     where  C  is  the  electro- 
static capacity  of  a  metal  disk  of  the  same   shape   and  size  as  the   aperture. 

The  variational  principle  for  ^,   was  successfully  applied  by  Levine   and 

r33l 
Schwinger  >•     ■■   for  the  diffraction  of  a  normally  incident  plane  wave  through 

a  circular  aperture.     The  correct  aperture  distribution  is  assumed  in  the 
formM 

(5.25)  9»i  =Z     a^(l-pV)""^/^  , 

n=l 

where  the  coefficients  a  have  yet  to  be  determined.  The  denominator  of 
equation  (9), 

(5.26)  F  =y^[k^52j^(p)52i^,(p')-V9l^(p)V'^^,(p')]G(p,p')dSdS', 
can  be  written  in  the  form 

1     /  /  r^"      f2  3^   (p)3!2     (p-)  -1 

(5.27)  F  =  i^/pdpyp'dp'/dY'[k''^^(p)i2^,(p')-  -3^ 2_ cos(Y-Y')J 


e 

X 


ik  \/p^-2pp'cos(T-?')+P'^ 


V^p^-2pp'cos(T-Y')+P'^ 
where  we  used  the  relationships  Vf5  (p)»V'^   ,(p')=  ^3^     -.  Y     cos  CY-Y')   and 


"^n       '^n' 


:V^^-. -- ' 


ikVP  -2pp'cos(Y-?')+p'' 
G       ® 


l/p^-2pp'cos(Y-Y')+p'^ 

Now  G  can  be  written  in  the  form 

ik\/p2-2pp'cos(Y-Y')+P'^  -,        00  ^ . 

(5.28)  ^  =  ^/   -^^  J^(XVp  -2pp'cos(Y-Y')+P''^). 

2nVp^-2pp'cos(Y-Y')+p«^  o    \/\  -k 
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Mien  we  substitute  equation  (28)  into  eqiation  (27)  and  use  the  addition  theorem 
for  Bessel  functions,  the- denominator  becomes 

(5.29)  F-  2it/°°-^^^/pdp  /p'dp'[kV(p)g    ,(p')J^(\p)J„(Xp') 

il-^J-^(Xp)J^(Xp.)J. 


3p        3p' 
__  _  i  .T   (\n  ) 

ap 


However,  the  expressiony  pdp  ^—  J,  (Xp )  when  integrated  by  parts  gives  us 


/pdp  ^  J,(xp)  =  P^^V^P)  ^  -y  ^n  apt  P  '^i^^'^^i  ^P 

o 

o 
a 


From  this  vie   see  that 

a         a  3^     352)    ,  _  a         a 

/pdp/p'dp'    gr^  3^  J-l(Xp)J^(XP')=  Xy  pdp/p'dp';?Jp)!2^/p')J^(Xp)J^(Xp«). 


Hence  equation  (29)  becomes 


(5.30)     F  «  2i 


°°X(k2-X^)  ^         ^ 


n/    ^^^1=^  dX/pdp/p'dp'J^(Xp)J^(Xp')9'j^(p)i3^,(p'), 
o    v^  -k  o         0 

yhich  can  be  written  finally  as 

00    r-s — p         a         a 
(5.31)     F  =  -2n/x\/x  -k^dxypdp/p'dp'52  (p)0   ,(p')Jo(Xp)Jo(Xp'). 

O  O         "Xl 

Equation  (31)   is  the  same  as  that  obtained  by  Levine  and  Schwinger.     However, 
they  used  a  non-integrable  kernel  in  the  integral  equation,   and  hence  the  valid- 
ity of  their  procedure  is   somewhat  doubtful.     Although  the  present  method  is  more 
difficult  than  Levine  and  Schwinger 's,   it  is  a  valid  procedure. 

Now,   if  the   series  expansion  (25),  where  the  individual  terms  are   of  the  form 
(1  -  p2/a2)^-l/2. 
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is  substituted  in  (31),  we  get  a  double  sum  of  terms; 

^5-32)       F^  =  -2n/'\VK?  dX   |  / pdp(l-pW)""^^^J  (Xp)| 

o  L  o  -^ 

[/pdp(l-p2/a2)'"-l/2j^(Xp)]      . 


o 

However 

a 


o     o  r.  MO                   o    y2                               -             a^J     T  ,^(Xa)2""-'-/^JTnH-l/2) 
ypdp(l-pVa'^)""-^/'^J^(Xp)=  a'^y  J^(Xasine)sinOcos^"»dO ElliZL. 


(Xa)"*^/2 


Substituting  this  into  equation  (32)  we  have 

(5.33)     F^.-„2"™((„4)[Tm4)a3/^HiiZ2!^!^Jli/£!l     V*?:?  dX, 

(Xa)"™ 

and  then  writing  X  =  kv  we  obtain  finally 
(5.3U)   F  -  -2Tta(2/ka)"''™"^  f?n+l/2)  [?m+l/2) 


nm 


V"(v'-l)'^'v-("«)j_,^^^3(k.v)J^^^^,(kav)d.. 


o        P    n    1  /P 

Moreover,  using  the  equation  ^  (p)   =  (l-p  /a  )  ~  '    ,  we  have 

-,  ft 

y  i2)^(p)dS  =  2nypdp(l-pW)""^/^ 


„      //^.    ^       2n^^„            2    r(l)frn+l/2) 
2tt  /     sinOcos     Od©  =   na      ^        '  ' — - 

o  F(n+3/2) 


2 

na 


> 


and  we  can  now  write  A,  (the  amplitude  in  the  forward  direction)  as 

r  °°       ^    12 
(5.35)  ^--^^ , 


2a       CO      00 
m=l  n=l 


^5*0     a  a 
*— -  ^-— _   mn  m  n 


where  the  coefficients  c       are  defined  by 

mn 
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(5.36)  c^  =  (2/ka)'"''''  Hin+1/2)  |Tn+l/2)f^(ka) 
and 

(5.37)  f^(ka)-/V.l)VV('"-^J^^^/2(kav)J^^^/2(^-)^-- 

Let  us  nov;  apply  small  variations  5a  to  the  true  values  of  a  .  In  view 
'^^  "  n       n 

of  the  stationary  character  of  (35)  we  then  find 
^    00   a     -A^   00 

^sn  Z  2HTI  =  2r  r  V^n  (-i>2,...). 

n=l  n=l 


On  the  other  hand,  it  follows  from  equation  (5)  that 

00   a^ 

z 

n=l 


(5.38)        A^^-ika^X  ^. 


Elimination  of  A,  thus  gives  the  following  infinite  system  of  linear  equations 


for  the  unknown  a  ' : 
n 


00 


(5.39)        Z.     c  a  -  i (m  =  l,2,...). 

^1  "^  "   ika(m+l/2) 

An  approximate  solution  of  (39)  was  obtained  by  Levine  and  Schwinger  by  assuming 
a  *  0  if  n  >  N  and  solving  a-,»*»,aj^,  from  the  first  N  equations  of  (39).  The 
corresponding  approximate  value  of  the  transmission  coefficient  becomes 

(5.U0,      4«.Ke  r   ^  a<«. 

n=l 

It  seems  worth  noting  that  the  integrals  (37)  can  be  expressed  in  terms  of 
Fjj  =  '^n'*'^''n  ^^®^®  Hn  ^^  Watson's  notation  for  the  Etruve  function.  The  symmetry 
between  the  real  and  imaginary  parts  of  f   were  not  recognized  by  Levine  and 
Schwinger,  although  their  expressions  are  easily  transformed  into  the  symmetrical 
form  by  a  partial  integration. 

T  Levine  and  Schwinger  apparently  overlooked  the  fact  that  their  coefficients  A  and 

[351  " 

D  are  simply  related  by  A  =  ikD  .  The  factor  C  of  Magnus*-  -■  is  thus  equal 
n^"'         •'nn  o^  ^ 

to  ik,  so  that  his  Table  I  provides  at  once  the  first  terms  of  the  power  series 

for  A  , 

n* 
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One  has 

ka     o         oa 
2  ?n 

oa     oa     ka.     o  16a  oa 

(5/Ul)         .3        o     I  c      T       o     I  c    2a 

oUa  2a    16a   o 

32a      Ua  Ua 

The  infinite  system  of  linear  equations  (39)  was  thoroughly  studied  by 
Magnus  '-    -^ .     He  showed  that,  for  sufficiently  small  values  of  ka,  the  solution 
a  is  unique  and  can  be  expanded  in  a  convergent  power  series  in  ka.  It  was 
also  shown  that  in  the  N —  approximation  the  exact  values  of  the  first  N  coeffi- 
cients of  the  power  series  for  a,,  ..,,a-,  are  obtained.  Explicit  recurrence 
formulae  were  given.  Special  attention  was  paid  to  the  limiting  form  of  (39) 
for  ka  ->oo.  Owing  to  an  error  of  sign  in  the  definition  of  the  imaginary 

part  of  c   ,  the  conjugate  values  should  be  taken  in  Magnus's  table.  The 

"^  r37l 

same  error  (and  others)  occurred  in  the  paper  by  Sommerfeld'*  •',  who  took 

2     2  1/2 

the  wrong  definition  of   (X  -k  )   '      when  0  <  X  <  k.     The  author  checked  Magnus's 

table  of  coefficients  and  found  complete  agreement  with  his  own  calculations, 
ej>:cept  for  the  last  column.  The  following  values  are  taken  from  the  author's 
paper  l-     -l;   i^re   give  one  term  in  addition  to  Magnus: 

^  -  -  ^  &  *  |(,a)^  .  |(ka,3.  ^i,.)K  ^(.a)5 

o±n 

i(ka)7  r      2  2x1 
For  the  aperture  distribution,  an  alternative  expansion,  which  is  somewhat 
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simpler^     than  the  Le vine -Schwinger-Sommerf eld  form  (39),  is  in  terns  of 
Legendre  functions.  Viz., 


n=o 


The  diffracted  amplitude  in  the  forward  direction  is  then  determined  by  the 

first  coefficient  alone: 

(5.'4U)  A^  =  _  I  ika^b^. 

Arplication  of  the  variational  principle  leads  to  the  following  infinite  system 

for  the  unknowns  b  : 
n 


(5.ii5) 
where 

(5.U6) 


00 


^^  ^mn\  =  iki  ^o  (^00  -  1'  ^mo-  °  i^  ™  >  o). 


n=o 


d       =   (6/ka)2    r(inO/2)|Tn>3/2) 
mn       '    f      '         ml  nl  ^mn^   ^ * 


g^(a)   =/  (v'-l)'/V2j^^^3/^(av)J2^^3/^(av)dv. 

In  this  case  the  various  approximations  to  the  transmission  coefficients  are 
given  by  the  ratio  of  two  determinants,   viz., 


fc\'?\   4-(N+l)     Uka     ^ 

Kt^.Ul)   t^         =  ^^     Imagm,  part  of 


=11 


•   •    •  g 


IN 


'Nl        •    •    •   %N 


in  ^ich  f or  N  =  0  the  upper  determinant  should  be  interpreted  as  unity.     It 
may  be  verified  that  equation  (U7)   gives  exactly  the  same   approximations  as 
equation  (39).     The  advantage  of  (U?)   is  in  its  explicit  analytical  form,  which 
invites  a  detailed  study  along  the  lines  of  Magnus's  paper. 

*  St=?rS^K  ^^^^l  ^^^  the  l^gendre-function  expansion  of  Levine  and  Schwin^er    r39l 
obtained  by  direct  integration  of  the   differential-integral  equation   (5.U).   "        '•  ^ 
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Only  slight  changes  in  the  preceding  analysis  are  necessary  in  order  to 
cover  the  second  boundary- value  problem.  First,  the  series  analogous  to  that 
of  (25)  represents  the  aperture  values  of  89(p/3z}  in  this  series  we  include 
a  term  with  n  -  0,  Secondly,  -^-k^/Za.   is  given  by  equation  (26)  modified  to 
include  the  terms  m,n  =  0.  Thirdly,  f   (a)  should  be  replaced  by 

(5.U8)        h^(a)  .  /u'-l)-'/'v-(-") J^^^/^(av)J^^_^/2(«-^^-' 

which  function  is  related  to  f       by   >-    •' 

ran 

(5.U9)  h     (a)   =   (m+n-2)f       -  a^f     (a). 

^  mn  mn  da     mn 


For  example. 


h     (a)   -^/Fn(t)dt, 


o 


(5.50)  hQ^(a)   =  I  .  ^/  ?^(t)dt  -  ^  Fj2a), 

Ua     o 

Ua     o  oa 

The  corresponding  first-order  approximation  of  the  transmission  coefficient, 
tp      ,  was  calculated  by  Miles'-     ■'    although   this  was  obtained  by  a  less  powerful 
variational  principle.     Miles  introduced  an  impedance  parameter  Z  =  R-iX,   and 
the   admittance  Y  =  Z~     =  G-iB,  vdiich  were  evaluated  for  constant  and  static- 
field  aperture  values  of  3{2„/3z  and  compared  with  the  rigorous  and  Kirchhoff 
values  of  the  transmission  coefficient  (tp  =  ReZ).     His  curves  are  represented 

by 

o 
where 
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(P,Q)  =  (1,0),  (0,1)  or  (1,2/kaii). 
Again,  the  simplest  analytical  approximations  of  the  transmission  coeffi- 
cient t^  (circular  aperture,  plane  wave  at  normal  incidence)  are  obtained  if 
we  start  from  an  expansion  in  I^gendre  functions.  We  now  assume  that  in  the 
aperture 


(5.51)       ^  =  7"  B 


9^9    °°      p.  (^C7/7) 


2   ^  .  '2n 


"=°  "   vCpVa' 


2.  2 


Then  by  equation  (21),   the   scattered  amplitude  becomes  A„  »  -a^B   ,  while 

c  o 

insertion  of  (51)  in  (22)  gives 

2 

(5.52)       ^=  2 

2a    00   00  » 


t— ^  ^^—^     mn  m  n 
m=o  n=o 


where 


Application  of  the  variational  principle  for  p     gives  the  following  infinite  system 


of  equations  for  the  unknowns  B  : 

n 


(5.5U)        21  D^3^  =  -|6^^  (m=  0,1,2,...). 
n=o 

The  successive  approximations  of  the  transmission  coefficient  are  then  simply 
given  in  close  analogy  with  equation  (U7),  by 
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11 


•  •  ^IN 


°N1   ••  •  °NN 


(5.55)      t 


(N+1)^  _  A. 
nka 


Ira 


G««  •  •  •  G-M 

00  ON 


No  •  •  •  ^NN 


It  should  be  noted  that  in  the  integrals  (37),  (U6),  (U8),  and  (53)  the  roots 

are  understood  in  the  sense 

(v2.1)i/2.  .i(l.v2)V2,  (v2.i)-V2  .  ,i(i.v2)-l/2  ^en  0  <  V  <  1. 


All  these  integrals  can  be  expressed  in  terms  of  f   and,  therefore,  in  terms 

mn  ' 


of  F  and  the  indefinite  integral  of  F  , 
n  o 
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71.     "Rigorous  Form  of  Bablnet's  Principle  in  Electromagnetic  Diffractlozi 

Theory, 

On  several  occasions  we  have  disctiseed  Bahinet's  principle  in  one  form 
or  another.   Only  recently  has  it  heen  possible  to  extend  this  principle  so 
as  to  "he  applicahle  to  rigorous  electromagnetic  diffraction  theory.  Oar 
representation  is  essentially  that  of  Cops  on  tj'*"'-^       .   In  what  follows, 
the  time  factor  exploit)  is  omitted;  as  "before,  k  denotes  the  wave  num"ber. 

Let  (f ,g)  denote  any  ar"bitrary  incident  field,  where  f  stands  for  the 
electric  field  vector  and  g  for  the  magnetic  field  vector.   It  is  assumed, 
therefore,  that  f  and  "g  satisfy  Jfexwell's  equations.   later  on  we  shall  use 
the  term  "complementary*'  incident  field.  This  is  the  field  defined  "by 
(^,f),  in  the  order  (electric,  magnetic)  vector.  As  is  well  known,  this 
complementary  field  also  satisfies  Maxwell's  equations.  The  complementary 
field  of  the  complementary  field  is  Identical  with  the  original  field  ex- 
cept for  sign. 

First  of  all  we  consider  the  diffraction  of  the  field  (f,g)  "by  a  per- 
fectly conducting  plane  screen  (finite  or  infinite)  of  zero  thickness. 
Secondly,  we  consider  the  diffraction  of  the  complementary  field  (-^,f ) 
"by  an  aperture  in  a  perfectly  conducting  plane  screen  of  zero  thickness j 
the  aperture  in  the  second  problem  is  of  the  same  size  and  shape  as 
the  screen  in  the  first  problem.  For  simplicity  we  call  these  two  diffraction 
problems  complementary  diffraction  problems.  The  rigorous  fora  of  Bcbinet's 
principle  asserts  that  the  solution  of  one  of  these  apparently  different 
problems  gives,  at  once,  the  solution  of  the  remaining  problem.  We  now 
turn  to  the  proof  of  this  statement,  and  to  Its  precise  form. 

In  the  first  problem  the  total  field  everywhere  in  spnce  is  given  by 
(f  +  E  ,  £•  +  H  ),  where  the  scattered  field  (E   .IT)  can  be  derived  from  the 

vector  potential  A^  of  the  currents  induced  in  the  screen  by  the  incident 

flow.  Let  I  denote  the  surface  current  density  vector.  Then 


-uo- 

H®  =  curl  a".    -iicE®  =  kV  +  grad  dir  A^. 

-A 
The  unknown  two-coraponent  vector  I,  defined  only  on  the  screen,  has 

to  satisfy  certain  integro-differential  equations  in  order  that  the  well- 
known  "boundary  conditions  shall  be  satisfied  at  the  surface  of  the  screen. 
The  superscript  s  will  he  omitted  for  all  quantities  evaltiated  at  the 
screen's  surface.  7or  example,  the  normal  component  of  the  total  magnetic 
field  must  vanish  at  the  screen.  This  requires  that  (  z  =  0  in  the  plane 
of  the  screen) 

(6  2)         ^      ^x 

r-^  -  r —  =  -  g  for  all  points  (x,y,o)  on  the  screen. 
ox  dy      z 

Similarly  the  x-component  of  the  total  electric  field  must  vanish  at  the 

screen.  Thus 


ikf^(i,y,o)  =  k\  +  |- 
X  X   ox 


3A    aA 
9x    3y 


3^A     ,   raA 

ox 


ay  Lay  "  ^zJ 


=  k^A  +  Aa  -  ^ 
X  X        oy 

=  k\  ^  AA^  *  i'rf:,  -  e/  (-.y.o), 

where  A  is  the  two-dimensional  laplace  operator  and  where  we  have  used 

the  facts  that  A  has  a  7ero  z-component  and  (f,?)  is  a  solution  of  ^fexwell•3 

equations.  The  last  equation  can  he  simplified  to 

2      A     ^^v 
(6.?)         k  A  +  Aa  =  r— ^    for  all  points  (x,y,o)  on  the  screen. 
X       X   oz 

Also,  from  the  condition  that  the  y-component  of  the  total  electric  field 
vanish  at  the  screen,  we  have 


-i.1  - 

2  ^^x 

(6.^)  k  A    +   AA     =  -  r —         for  all  points    (x,y,o)  on  the  screen, 

y  y  oz 

FinP-lly,  'bearing  in  mind  that     on  the  screen     we  have 

/.         ilcr  /         ikr 

(6.0  A     =i/l     2 ds  A     =i    /l     2 dS, 

X      cy     X     r  •  y      c  y     y     r  » 

where   r^  =  (x-x«)^  +  (y-y')^,  dS  =  dx'dy',  "l  =  "iCx'.y"),  we  see  that 
eanptions  (2)  throti^h  (5)  constitute  a  set  of  inte^ro-differential  relations 
for  the  ijnVnown  cnrrent  density  I. 

By  physical  intuition  we  expect  these  relations  to  have  at  least  one 
"admissihle"  solution  I  =  (I  ,  I  )  satisfying  all  physical  req^uirements  as 
to  singularities  possibly  occTirring  at  the  ed^e  of  the  screen.   It  is  not 
known  whether  the  assumption  of  absolute  integrability  of  I  over  the  screen 
would  entail  a  unique  admissible  solution.  The  integrals  in  (5)  cannot  be 
proper  Riemann  integrals;  they  are  improper  because  of  sing\ilarities  of  I 
and  I  at  the  edge.  Recent  work  of  Meixner*-  •" ,  Maue"-  •* ,  and  others  makes  it 

plausiblp  th^^t  the  component  of  I  tangential  to  the  ed^e  becomes  infinitely 

-1  ^2  ->■  1/2 

^pTf(^   a«<  D     r>n(1  that  the  component  of  I  normal  to  the  edge  vanishes  as  D  '  , 

where  D  is  the  distance  to  the  edge.  Similar  properties  hold  for  the  field 
vectors  themselves,  although  it  is  not  clear  at  present  what  conditions  are 
necessnry  and/or  sufficient  for  a  unique  physically  acceptable  solution. 

Ve  now  consider  the  complementary  diffraction  problem.  There  is  need 
for  a  proper  distinction  >ietween  the  fields  in  front  of  and  behind  the  aper- 
ture. Let  (E  ,  F  )  denote  the  total  field  in  the  illuminated  space  If  there 
o   o 


is  no  hole  in  the  screen.  I'or  exf>mple 

y.-z)  1 

z  ^  0. 


5!^^  =  -gj^(x.y.z)  +  gy(3f,y,-z) 


V  "  ^y^^'J^'^)  •"  fy(^.7.-z) 
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I«t  (E^,  H^)  denote  the  total  field  for  z  ^  0  in  the  presence  of  the  hole, 
and  let  (B^,  Hg)  denote  the  total  field  for  z  >  0  (that  is,  "behind  the 

aperture).  Further  let  K  be  any  two-component  vector  defined  in  the  open- 

^d 
ing.  Befine  the  vector  T     hy 

{6,6)  y*  =  -i/^K2 — ds. 

c  V    r 

Then  we  will  show  that  when  K  =  (Z  ,  K  )  satisfies  the  proper  conditions, 

the  fields  wn  he  expressed  as  follows: 

/^       -»        -»  -'d 

'         E,   =  E     -  curl  F  /         /     ^  «\ 

1         0  ^  [        (z  <  0) 

-IkST   =  -ikH*   +  k  F*^  +  prad  div  f' 
F     =  curl  r" 

P-,        ^.      ^   (^>o). 

+  IICH2  =  }r:'^     +  grsd  dlv  F^ 

First  of  all,  so  lon^  as  r  is  Inte^rshle,  the  field  defined  "by  the  preceding 

-•    -»d 
equations  In  terms  of  K  vie  F  satisfies  Maxwell's  equations  in  the  half- 

spacRs  z  <r'  0  and  z  >  0,  and  it  satisfies  the  appropriate  "boundary  conditions 

at  the  screen  (i.e.,  thft  the  tangential  electric  and  normal  magnetic  fields 

vanish).  Furthermore,  no  matter  what  K  ig,  the  tangential  electric  and  normal 

magnetic  fields  are  automatically  continuous  in  the  aperture. 

The  only  conditions  that  are  not  automatically  satisfied  are  that  the 
normal  E  component  and  the  tangential  H  component  be  continuous  in  the  aper- 
ture.  In  fact  these  components  are  easily  seen  to  be  identical  with  the 
corresponding  components  of  the  undisturbed  incident  field.  Again  dropping 
the  superscript  d  when  we  refer  to  values  on  the  apertiire,  we  find  that  these 
conditions  are 

aF   aF 

(6.8)  -r-^  -  -r —  =  -g    ,      for  all  points    (x,y,o)    on   the   opening. 

ox         oy  z 


-1*3- 


(6.9)  ^\-^^^.  =  ai^^ 


for  all  points    (x,y,o)  on  the 


(6.10)  k%  *  AJy  =  -  ^   I  ''^''"^- 

It   is  remprVs"Me  that  these  eqtifitlons  are  precisely  the  same  as  those  for 
A   ,  A     In  the  first  pro'blem.     AssTimin^  that  either  system  has  a  \inique  solu- 

tion,  we  see  that  the  two  mutijally  complementary  pro^blems  are  simply  Identical i 

— *■         J-       — i  _i  -»-£[         -^  g 

™  =  A,  K  =  -I;  hence  T     =  A  .  Moreover,  If  we  introduce  the  notation 
(^2*  ^9^  -    ^^   •  ^  ^'^^   *^®  diffracted  field  behind  the  aperture,  we  have 

(6.11)  B'^  =  H^  and  H*^  =  ^*        (z  >  O)^ 

and  this  is  '^'binet's  principle  in  its  rlp-orous  form  for  diffraction  "by 
plane  perfectly  conducting  screens  and  apertures. 

By  analogy  to  the  first  prohlem,  K   is  called  the  magnetic  current 
density.   Its  components  show  the  same  hehavior  at  the  edge  as  does  the 

electric  current  density  In  the  complementary  pro'blem.  The  vector  T     is 

called  the  magnetic  vector  DOtentlal. 

In  concluding  this  section.  It  nay  he  noted  that  taking  the  comple- 
mentary field  [transformation  (f,g)  ->  (-^t^B  of  a  plane  wave  is  equi- 
valent to  rotating  the  plane  of  polarization  through  a  right  angle  counter- 
clockwise, looking  In  the  direction  of  propagation. 


-  kk  - 

VII.  Diffraction  by  a  Small  Aperttire  In  a  Perfectly  Conducting  Plane  Screen 

In  this  section  we  shall  comment  upon  Bethe's  '-  ^    theory  of  electro- 
magnetic diffraction  by  small  apertures.   It  appears  L  -1  that  Bethe's  first- 
order  approximation  Is  fundamentally  Incorrect  with  respect  to  the  field 
near  the  aperture. 

Let  E  ,  H  denote  the  Incident  field,  and  let  1  ,  H  denote  the  field 

In  the  illuminated  half-space  (z  <  O)  if  thpre  is  no  eperture  in  the  screen 

(z  =  O).  Then 

E   =  E^C?)  -  E^(-z) 
ox    X       X 

E   =  E*(z)  -  E*(-z) 
oy   y     y 

E   =  E^(:^)  +  E^(-z) 

0  3      Z         Z 


ox      X  JL 

H   =  H^(z)  +  H-t(-z) 
oy   y     f 

H   =  H*(z)  -  H^(-z) 

02      Z         Z      ) 

where  z  <  0  and  Where  we  have  omitted  explicit  reference  to  the  x,y-coordin- 
ates. 

As  we  have  seen  in  the  preceding  section,  the  diffracted  field  can  be 
derived  from  fictitious  magnetic  currents  (and  charges)  in  the  aperture: 

"H,  =^  +  IkE  +  (l/k)  grad  dlv  E  (in  front  of  aperture) 

1     0 

(7.1)  I     _ 

^2  =  ^^'•^  ^  (2  >  0) 

-»     -^   ,  ,  V         :^  (behind  aperture) 

H2  =  -IkT  -  (i/k)  grr.d  dlv  F  ^  > 

where  7   la  the  magnetic  vector  potential  given  in  terms  of  the  currents  K 
by  means  of 


-i/ 


-  /i5  - 
Ikr 

(7.2)    ^  =  -^y       kv*s- 

apertvire 

The  magnetic  chprpe  density,  y   ,    is  fo"und  from 
(^.^)  div^  =  _^    ^.     --Z    .  i^c^^ 

and  the  scalPr  magnetic  potential    w  satisfies  the  equation 

/  ikr  _^ 

(-7.1)  ^  =  y  r7^dS  =  idiTF. 

flpertijre 
In  order  that   these  formulae  hold  it   is  necessary  that  the  component 

of  K  normal  to  the  ed^e  of   the  screen  vanlBhes  at  the  edgeC   cf .  Boiwkamp '-  ■'•), 

For  arbitrary  K  the  tangential  electric  and  normal  magnetic  components 

are  automstically  continuous   in  the  aperture.     We  hsve  for  these  components 

('.5)  t=|;GSx^:.  y.j^iX^)^ 

where  script  letters  refer  to  valines  in  the  aperture,  and  n  is  the  unit 
vector  in  the  positive  z -direction. 

Conditions  for  K  are  ohtained  by  requiring  that  the  normal  electric 
and  tangential  magnetic  comrjonents  are  continuous  In  the  aperture.  Thus  we  have 


£^  =  lim  E^^  =  lim  E^^  =  p 

3t     -   llm  H,  =  lim  F   «  r  H    (and  same  for  y-comp.). 
X       ix       23t   2  ox 


In  the  opening,  therefore. 


2^  0 

Consequently,  as  a  by-product  we  get  the  theorem:  in  the  aperture  the  values 
of  the  tangential  magnetic  and  normal  electric  field  components  are  exactly 
equal  to  the  values  of  the  corresDonding  components  of  the  undisttwrbed  inci- 
dent field. 


-  liS  - 


As  we  have  seen  In  the  preceding  section,    the  ahove  requirements   lead 
to   the  system  of    Integra l-differemtial  equntions 

ilcr 


/ 


(-^.7)     / 


X  c  /     X     r 

,     r  ikr 

<J    =  .  i  /  K     2_- 
-'7  c  /      y      r 


dS 
1 


X  X  9z 

^y  y  3z 


> 


^    ax  ay         *^z 


T=T(x.y) 

"?=  K(x',y') 
dS  =  dx'dy' 

r^-  (x-x')^+  (y-y')^ 

(x,y,o)  any  point  in  aperture 
(x'.y'.o)  same 


So  fpr  our  fomrilPtlon  is  general.  We  next  consider  the  case  of  long 
waves  (k  ->  o),  and  assume  thet  K   can  "be  expanded  in  a  power  series  of  k. 
I«t  K     pnd  K     denote  the  psrtB  of  K  of  relative  order  zero  and  one  respect- 
ively.  It  enpears  that  relative  order  is  the  same  as  ahsolute  order.  Thus 
assume 

"k  ""k®  ■*-'V-  +   O(k^). 
where  K     =  0(l),  &  «  0(k).  Then,  expanding  the  exponential  function  in 
powers  of  k,  we  hprrf^ 

The  third  integral,  of  order  k,  apparently  is  independent  of  x  and  y.  The  re- 
sult of  differentiation  of  this  term  with  respect  to  x  or  y  is  identically  zero, 
We  need  not  retain  this  term  in  the  expressions  for  the  fields  in  otit  order 
of  approximation.  Therefore  J  =  T  +  J   +  const. k  +  0(k  ),  where 


(7.8b) 


^°._i/2_^S  (order  1) 

J^S-i/^dS.  (order  k) 


with  the  same  convention  ps  to  the  suiDerscripts   7ero  and  one  as  "before. 


-  ij?  _ 


We  have  C^ ^  A'j®  +  A!f  +  0(k  ).  Therefore, except  for  terms  of  order  k  , 
we  ina7  replace  k  v  +AT  by  A^  .  Eqtistlons  (?)  then  reduce  to  (8a)  and  the 
following 


^^. 


(7.8T>)S^^  y 


9zJ 


k=o 


*~  k  =0 


[•&b  -<-fa-%i. 


'k=o 


k=o. 


Let  us  evfllupte  the  rif-ht^hand  sides  of  equations    (8b)  for  a  plane- 
polFTlrer"    o'Kl-ianply  inelrient  electromn^netlc  wave.  We  choose  the  xz -plane 

es  plane  of   incidence  and  call  0     the  angle  "between  the  z-axis  and  direction 
of  incidence.     The  phase  function  of  the  incident  wave  therefore  is 
exp[ik(xsinO  +zco86    3*     ^'urther,    let  f)     denote   the  angle  between  E     and  the 
xz-plane.      Then,   In  en  ohvious   component  notation, 

E     =    (co9^  C08«   ,sln^   ,-co80  sine   )exp[lk(xsin©  +2Cos©   T] 

H     =    (-8ln0  cosG   ,cos«5    ,8in(5  sinO   )exp[ik(xsin©  +zco3©   ff  , 

O  0  0  0  0  *^  0  0 

Accordingly,  in  this  case  we  have 

-  r-^  =  -I'csiniJ  cose  expdkxsin©  ) 
oz  0    0  0 

- — m     IVcos^  cos  ©  expdkxsln©  ) 

07,  0  0  O 

£,     =  -co3(i  sln©^  expd'ocsln©   ), 


and  eou?>tlon8    (8h)  'heoome 


-  1j8  - 


(•7.9)  \     A^°  =  0  A'f^   =  icos0^co8^e^ 

Note  that  if  we  had  replaced  -  — L     ,  +  r —  ,  <^  in  equetions  (?) 

doz  z 

z 

by  their  corresponding  constant  values  at  the  origin  of  coordinates (which 
may  be   the  center  of  a  circular  aperture,   for  example)  we  would  have  obtained 
the  apTiroximfite  equPtions 

A^    =  -ik8in{5  cosO 

^X  ^0  0 

(7.10)      /  A^    =  iicco80  cos^e 

»  -^  y  0  0 

r-^  -  T —  =  -cos^  sine  , 

3x         3y  0         0  ' 

and  consequently,  bearing  in  mind  that     J  =    J     ■*■   )     ,  we  would  have  obtained 
eqs.    (9)  except  -for  the  Ifist  one: 

r— ^  -  -r—  =  0  instead  of  =  -ixcosjj  sin  ©  . 
ax  dy  0  o 

Faufftlons    (10)  are  essentially  those  of  Bethe^  -^  and  Copson^-^.     They  can, 
therefore,   only  le^d  to  a  correct  '^^°  term;   the   term  "^      is  necessarily  wrong 
in  their  approximation. 

It  is  possible  to  eliminate  derivatives  with  respect  to  k.   by  introducing 
derivatives  with  respect  to  x,y: 


-co8|j  Sine   [i+nacsinell««5  (0,0, 0)  +  x  r—  (0J0,0). 
o         0  cT        z  ox 


We  then  have  expressed  all  constants   in  terms   of  center  values.     More  generally, 
we  assume  ^  4 

6^(^.y. 0)^5^(0,0,0)  +   X  -gfto.0.0)  +  y  ^(0.0.0), 

Our  final  equations   thus  become,   for  any  plane-wave  excitation. 
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(7.11)  A'J°  =  0-         (7.1^)  A^^  =  -P  •. 

(7.12)  A^Jy  =  oj       (7.15)  A'fy  =  q  i 

as*    as*    aE*    as* 

where  P  =  r;;^^  .  q  =  --^,  R  =  _-Z.g^^y„j;*, 

All  these  eqw-tlone  are  evaluated  et  the  center  of  coordinates.  The  constants 
h&To  the  following  orders  of  magnitude: 

^   =  0(1);  P.Q.E.S.  =  0(k). 
This  theory  will  now  he  applied  to  the  case  of  a  circular  aperture  of 
smPll  radius  a  (l^  «  i).  p,,  ,,,,  ^^^^^  ^^  ^^^  ^^^^  complicated  integrals 
of  the  type 


If 


')' 


f  =  1  ,  then  G  =  1  and  AG  =  0 

f  =  X '  ,  then  G  =  -  X  and  A  6  =  0 

f  =  x'y'  ,  then  G  =  ^  xy  and  AG  =  0 

f  =  x'^  ,  then  G  =  ^(^^+5x^-y^)  and  ^G  =  |  . 


(l^or  more  general  results,  see  Bouwkamp '-'*'^*^'^  .)  In  virtue  of  this  we  find  that 
1     P^-^^x'-^p^y'-hp^x'y'^p^Cx'^.y'^) 


c      X 


2\r2      |2      T 
T     ya  —I  •    -y ' 


1      0       Vqiy'-^V'-^q^3r'y'-t-q^(y'^-x'^) 

■•    ■"   K      s    ^^™™— — ^^  "  " "  "  ■■■  ^——^1  ■■».»■ 

c     y 


2\/2      ,2      J" 
T    ya  -x'   -y' 


solves    (ll)  and    (l2)    [with   10  arhitrary  parnmeters    (p  and  q2 •     Furthermore,    on 
evaluation  of  the  relevant   inte#:rals  we  find  that 


and  thiis 
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czo  1  1  3  3,22. 

^x  "  V  2  V  2  V  8  ^3^  *  8  ^^^^  -^  ^ 

^y  =  V  I  ^1=^  *  2  V  *  8  ^3^  *  8  ^2.(y^-^^) 


^-df'  2^VP2^  ""  8   C<l3+2p^)7  -   (p3+2p^)^. 
This  should  eqiial  V,   therefore  q,  »  -2p^f  P,  *  -2^2^,  Q-p-Po  ~  ^*     '''liea  *«  !>«▼• 

1  ^o    PQ-^Pi^^Pgy-^^^^^/j  <'^-y^) 


-  -  K 
(7.17) 


^l^rr^T^ 


1  o  ^  '^o*\y'^2^-^'^ii'^%  ^^^"^^ 
~  c   y  ~      oif  o   o   2  * 

7 


2\/  2     2     2 
Tf    ya  -X  -y 


whlc>i  contPlns  7  Indepfindent  pprpmeters  (8  pprameters  p,  q  with  one  relation, 
q  _p  =  2V,  -hfttween  them),  is  a  .general  solution  of  equations  (ll),(l2),  (13) 
end  the  first  equation  (8a).  Equation  (17)  Is  presuinably  the  general  solution 
with  the  requirement  that  K  "be  absolutely  integrable  over  the  aperture. 

It  is  clear  that  additional  conditions  are  necessary  in  order  to  find 
a  unique  solution. 

Of  course,  the  charge  density  should  he  integrahle  over  the  aperture. 
In  our  case,  this  mesns  that  xK  +  yIC°  should  contain  a  factor  a  -x  -y  , 

^         y 

liqulvalently,   the  rndial  component  K     of  K  must  vanish  at  p  =  a.      If  we 
assume 


(xKW)  %^-x^-j^     £    (d^+d^x^-d2y)(a2-iV), 


we  find  that  all  coefficients  (also  d  ,d-,d„)  are  to  be  zero,  except  for  the 

0    1    £ 

one  relation  V^'^o  ~   ^*  *•  ^'®  then  left  with  a  unique  solution 
ex    -.  r-T — -  *     c  y 


„2\/  2     2         °     ^       „2\r2~ 
TT    ya.  -p  n    ya  -^ 

or  In  vector  notation 


2 
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(7.18) 


2rr      ya  -p 


■NTote  that  dlv  K'  =0. 

In  the  same  way  we  determine  K  .  We  split  this  vector  into  four  simpler 

ones  which  can  he  more  simply  treated: 

First  Fart: 

In  what  follows  the  notation  (a  ->h)  means  "contrihution  to  a  is  h". 


ife(l)  _  2P(-2a^-4-2xW) 

""^       ~     ^Z\l  2     2     2       ' 
yt    ya  -I  -y 


c    y 


2P  xy 


3n2V?:;;v  * 


^^  ^    ^  (20a2.9x2-3y')  , 


^  ->  J  ^  ^ ; 


at^^ 

-  P 

Sx       " 

39^ 

ay 

0. 

Second  Part: 

1  k^2^=  . 

2Q  xy 

" 

C      X 

^^V.^- 

^'-7^ 

a7,  ^  0  i 


1  ^(2)  _  _  ;0(-2a^ti^<-2y^) 
o     y       " 


y?][J^?I? 


^.^   E»-^, 


<T,  ^  -  I:  (20.2.3,2-9/)  ; 


AT,^  0. 


A"Tv  ^   4  1 


ax  ay 
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Third  Part: 


1  ^('^)  =         ^5  xy 
c     X        ~ 


3" 


^yXX? 


^y^~     («a^+3x2-37^) 


A'Ty  ^  °   i 


y  X 


-^-R  X  . 


'"'ourth  Part; 

1  J^)       25  (^^-tV) 


%-^|sxy; 


It  can  "be  verified  that   the  corresponding  charge  density  contributions  are 
determined  by 

i  .._^(1) 2Px 


dlv  k' 


.'\/X, 


1  4i,t'2'  -  -  -^ 


v']/^, 


1  ^4    -r(3)  2Ry 

—  dlv  k         =  ^ 

c 


2  \l  ?     2 

TT     ya  -p 


c 


2\/  2     ? 

TT      V«    -p 


and  hence 
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(7.19) 


1  ,i,t^  „  ^   (P-s)x^(M)y 

c         z 


-    sin^ 


By  vlrtwe  of  Maxwell's  eq-uRtlons  for  the  incident  field,  we  have 

P  _  S  =r  -ik  H*  (0,0,0) 

R  _  Q,  =  >lk  H^  (0,0,0). 
7 

Moreover  div  Z**  =  0,  so  that  the  magnetic  charge  density  becomes 


-j.^          ^.'1(0,0.0) 
(7.20)    7  =  ^  div(K°+K^)  =  -  -^       °   


V^ 


(Bethe)  . 


li^l 


Bethe 's  expression  for  —  K  (hieK^^)  is  wholly  incorrect,  although  div  Z^ 


C  H 

is  correct.  According  to  Bethe 


\  =  %   Va^?  (0,0.0). 


"hut  the  correct  value  is 


\<-^\\^±2 


(P-2S)(a^-x^-2y^)-i-(2R-Q)«y 


>       \/ 2  2  2 
^    Va  -3C  -y 


-  Z^  =  (K„) 
c  y    H  y 


3" 


(P-2S  )xyf  (2R-Q.)  (a^-2x^-y^) 


v; 


2  2  2 
a  -X  -y 


Tho  correction  is  precisely  a  divergence-free  vector.  This  difference  can 
be  written  as  the  curl  of  a  certain  vector.   In  fact,  one  finds 


-  K"  =  K„  +  -2:r  curl  p  ^a^-x^-y^  { (P-2S )y  +  ((i-2R)xj  J 

o  r    si  -y      o   9  (   ^^^      SE^      dE^      aE^  -)-l 


3TT 


Using  vector  notation  we  have  the  following  expression,  which  can  he  applied 
to  plane-wave  excitation: 


>  5ii  - 

(7.21)   ^t^  =  -^  [Va%^{21^^+^x  gr8d(^  •  B^)} 

*-^  .  jikH^xn-grad(^.  B^)]  j^^  H  ^ 

fa  -p 

where  all  qiiantitles  refer  to  values  at  the  center  of  the  apertvire  and  the 
gradients  are  taken  in  the  tangential  direction  of  the  apertiore;  p  is  the 
radius  vector  from  the  center  of  the  circular  aperture, 

?!auation  (2l)  shows  that  a  vanishes  at  the  edge  aa   y*  "P  ^^^  thJBit 

0 

Kj   >»ecome8  Infinitely  large  as  1/  ]/&   -p^  s ince^n  X  p jhas  th6  direction  of  0. 
The  complete  expressions  for  our  plane-wave-incidence  case  hecome 


-  Ip-  =   ^ya^-p^  [lsin0  cos©  co8^co80  (l-  'Isln^O  )sin^'2 
(7.22) 

J  K^  =  ^  -p=  Csin)«^cose^(a2_  i  p2)sin|!^cos0 ^(a^-  |p2) 

+    (p^  -  ^%in^Q^  ]  COS03. 
^or  the  total  field  in  the  aperture  we  get 


5.  =  - 


2cos0  sin© 
0         o 


p  Tf  1/72     2 

Va  -p 

81k    ^^ r.._^   — n    ^«2     1   _2, 


-  -^    i7T=^  &in?5^cos©^(a  -  I  p  )sin0 
Ka  -p 

(-'.2")  +  cos«l^  {(a^-  I  p2)  +    (p2-  i  a2)8in^0^  cos^U   +  0(\c%.^) 


^tk  "  ^Va%^  G.sin0  cos©  co9^co80   (l-  ^in\)sinO  -f  o(k^a^) 


(j        -^ 


Jt^  =  P  [8ln<5^cos©^co8|!Ucos0^sinS  +  0(ka) 

-P 


.vr2~2  "  °    ° 

TT  Va  -j: 


where 


^.   f      ^  .      -^^    are  the  undlfltur>^ed  incident  field. 

I  p  p 
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For  the  case  of  normal  incidence  (Q  =0,  0  =0)  this  result  reduces  to 

0      0 

the  first  few  tern«  of  an  earlier  eaqpresslon    (BoviwkB.mp '-  -' ). 

Once  we  know  small-aperture  approximations  for  K,  (see  eqs.  (18)  and  (21)) 

we  can  try  to  ©Taluate  the  vector  potential  F  to  the  same  order  of  approxl- 

nation  hy   menns  of  (2)  and  ultimately  to  arrive  at  the  corresponding  field  hy 

means  of  (l).  This  process  is  easy  for  the  field  at  large  distances  from  the 

aperture  (in  the  wave  zone).  Applying  known  methods  Ue  find  the  asymptotic 

representations 

.3  ikr 

T^-v^  ~-  |P-S+  |ikV8in«8in0^7" 

--^Ca.|vslnesinO^ 

=  -^[>^|vsln«cosOV-  ' 
where  o,  =  H  (0,0,0),  P  =  H^(0,0,0),  and  P,  Q,  B,  S  have  the  same  meaning  es 

^         y 

before ,  while  r,  Q,  f>   are  spherical  coordinates. 
In  vector  notation,  one  has 

where  r  is  a  unit  vector  in  the  r-direction;  the  values  of  H  and  E  are  those 
o  0     0 

taken  at  the  center  of  the  aperture.     Using  CTirl(i  xr   )  =  0(7r),  we  find  that  in 

0   0       ' 

the  wave  zone  E  =  curl  F  +  0(p),  so  that 

2  3      ^  ikr 

(■^.25)    E,^  ^:^'r  T   r2H  +E  X?  1 

2    3TTo'-ooo-'r 

^,23  ^    ikr 

H-^  r  xE  ~  ~-  T  X   \%  +T   x2H  1- . 

2    0       3"o'-ooo-'r 

These  relations  ere  those  of  Bethe,  Although  his  expression  for  K  was  wrong, 

the  far-field  expressions  tutn  out  to  he  right.  The  explanation  is  that  the 

correction  term  in  ^  Is  a  divergence-free  vector. 
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The  field  (25)  'behind  the  eperture  In  the  wave  zone  Is  exactly  the  stun 

of  an  electric  dlpole  (moment  P  )  and  a  magnetic  dlpole  (moment  P  )  hoth  lo- 

6  in  ^ 

Gated  at   the  center  of  the  sperture: 

P     =  +  "TT  E  (normal  to  plane  of  the  aperture) 

(7.26) 

P     =  -  •^jfj"  ^  (in  the  plane  of  the  Bpertnre). 

Note  that  the  sign  of  P     Is   Incorrect  In  Bethe's  paper*-    -* . 

We  now  turn  to  the  transmission  cross-section  of  the  small  aperture. 
The  procedure  for  calculating  Is  well  known: 
!E|2  =  k2lr^:j|2  =  k2(j,2^  J.2J 

3        Ikr 
T     =  cose[j  co^F  sinS    =  -       Jl^ [ocos^PslnScos© 

J  J  =  -F  sin^r  cosjj  =  -     ^^  2^  [Iaain0+poo8^  |  Vslne] 

y%^|i|2  ^  I6a^    J^  Da^+P^)(l-^cos2oV.  I  V^sm^O" 
0  9Tr        r 

y     d^y     sln«!E|    d©  =  -^;;^ ^  L  3  («  +P  )  +  3  ^  J 

o  o 

Let  S         denote  the  time-average    (indicated  hy  a  har)  of  the  Poyntlng  energy 

To" 

flux;   then 


27TT 

This  agrees  with  Bethe's  expression    (we  have  corrected  a  trivial  misprint   in 
his  paper).      In  the  plane-wave  case  one  has 


AH^  +  e2  =  2!ii(8ln^0  cos^e  +cos^jj^)+cos^0  sin^«   H  . 


0  0 
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The  right-hand  side  hecomet 


8[l  +  I  sm^  ojj  for  0^=  0. 

8  co«    «  for  |!5^  =  X  . 

O  O  <£ 


Therefore 


„               8c       6.  ^       2- 
S  ,       = a  k  COB  Q 

27TT 


(".25) 

S 

2'?TT 


-  Be       6,  ^    n^  1  -J    2-    -| 


The  energy  flm  of  the  incident  wave  per  unit  of  cross-sectional  area  eqriAls 
c/grr.      Hence  for  the  cross-eectionSj 

cr.     =  — —  a  k  cos  Q 
-L         2^  o 


(7.''0) 


^      -^aVCl.l.lA     1. 


Tor  the  general  case  ^ 

(7.:^1)        cr    =  1^  Cl*   (Jcos\-8in2^^)sin\:. 

Tor  tmpolsrized  incident  vireves    (natural  "light    •), 

(7.32)  a"''-^     B-i=in\j. 

Equations  (28)  through  (32)  are  d^re  to  Eethe, 

A  few  words  must  he  said  ahout  the  field  near  the  aperture  (the  qyasi- 
static  field),  ^or  this  case  we  have  to  calculate  F  for  small  values  of  ka 
and  ke.  This  is  most  effectively  done  in  spheroidal  coordinates  u,  v  related 
to  the  cylindrical  coordinates  hy 

p     =  a    (l-u  )(l+v   ),   2  =  auv,   0^u<l,  -oo<v<   oo  • 
(See,  for  example,   Bottwkamp       "^  . )  Very  tedious  and  complicated  calculations 
are  required,    though.     Therefore  we  shall  give  final  results  only.     Note  the 
hrnnch  0  <  arc  cot  v  <  tt, 

Totpl  electric  field  near  aperture,   for  hoth  z  <  0  and  z  >  0: 
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^^     rQ(varccot  v-l)  +  ^^  ^  1 


gmc 


TTa(u2+v^)(l+v^) 


1     =  ^Si 
y         TT 


P(vprccot  v-l)  + 


2S-P 
3(l+v2) 


2T3.y 


Tia(u^+v^)(l+v^) 


2V 

1     =  — 

Z  TT 


r    .      V  1 

srccot  V  -  ~2 2 

u  +v    J 


+  -  (Ex+Sy)[?>rccot  v  -  r  - 


k^ 


1+v^   3(nV)(l+v^) 


^v(Q7+Py) 


3n(u2+v^)(l+v^) 


"For  the  totnl  mprr\ptic   field: 


IT  =  — '  jarccot  v  -   ^'J- 


H:c(ai»Py) 


^   TT •i^?-'   «.2/  2^2wt.  2x2   J 

1+v     Tta  (u  +▼  Hl+v  ; 


H 


y     TT 


f  E 


srccot  V  - 


rD- 


»vy(n3e-pyj 


2/  2.  2n/,^2n2  > 


l+v'^    Tia^(u^+T'')(l+v'') 


^n(axfpy) 


H  =  - 
^  TTa(n2+v2)(l+v2)  ■ 


Note  that  on  the  aperture  v  =  0,  u  =  yl-P  /*   ^>^^  on  the  screen  u  =  0, 
V  =  ±  Vp  /a  -1   (for  z  =  ±0).  By  calculating  the  Jump  of  H  across  u  =  0 
we  are  able  to  calculate  the  electric  current  density  in  the  screen  near  the 
edge.  This  is  left  to  the  reader.  For  the  case  of  normal  incidence,  see 
Bouwkamp  *. 
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VIII.     On  CoTJson's  Theory  of  Dlffrpction 

^^  ^ 

In  a  recent  jjpper.  Professor  Copson  '-  -•  hp-s  commented  upon  the  ■  note  "by   the 
reTlewer**  that  w?8  added  to  my  review  I-  -^  of  Copson 'o  efirller  paper  entitled 

"An  IntefTFl-eatiPtion  method  of  solving^  plane  diffraction  problems"  ^  -' ,  As  is 

C52I 
obvious  from  Copson's  onswer  as  well  as  from  two  papers  "by  Miles    ■*  ,  my  comments 

have  led  to  some  confusion  as  to  the  "boundary  conditions  to  "be  satisfied  on  the 

rim  of  "the  diffrr-ctinp  obstacle.   In  order  that  any  possible  misinterpretation  of 

my  criticism  mpy  be  avoided,  it  seems  worth  while  to  discuss  these  questions  in 

greater  detail  than  is  permitted  in  a  short  review.  This  also  presents  the 

opportunity  to  stress  how  carefully  Copson's  differential-integral  equations  must 

be  handled  in  practical  applications;  we  now  know  that  his  solutions  for  the  small 

circular  disk  and  apertrtre  are  in  error*. 

To  present  a  clear  picture  of  the  problem,,!  shrll  quote  the  essence  of  Copson's 

theorem  h   as  well  as  my  comments  published  in  Mathematical  Reviews.  The  theorem 

in  question  is  '-  J   • 


♦  It  should  be  noted  that  Copson  now  agrees  fully  with  the  theory  outlined  in 
this  section  '-  -^  . 

•*  For  the  discussion  in  this  chapter  we  have  preserved  the  origine.1  notation  of 
Copson.  Hla  f piTrilation  in  terms  of  scalar  functions  is  essentially  equivalent 
to  that  given  in  Section  VII  with  the  following  exceptions: 

(1)  In  Section  VII,  the  time  factor  was  taken  as  e~   ;  Copson  Mses   e    , 

(2)  In  Section  VII  the  exciting  wpvp  is  incident  from  the  left  while  here 
incidence  from  the  right  is  assumed. 

The  functions  v  and  u  are  the  x  and  y  rectangular  components  respectively  of 
the  vector  potential  7;  «  is  the  scalar  potential  \jf   , 

The  functions  e  and  -e  ere  proportional  to  the  x  pnd  y  rectangular  components 

y     ^ 

respectively  of  the  magnetic  current  density  K;  b  is  proportional  to  the  magnetic 
surface  charge  density  »[  • 

When  thesp  differences  are  taken  into  account  the  formulation  in  terms  of  vector 
functions  given  in  Section  VII  is  equivalent  to  that  given  here.  The  correspondence 
between  the  equations  may  be  seen  from  the  following  table, 

Equations  in  Sectien  VII  Corresponding  Eqxgtions  in  Section  VIII 
7,1  8.1  and  8.6 

7.2  and  7.1^  8,2  and  8,3 

7.3  8,^ 

7.7      .  8.5 
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Let  an  electromagnetic  field  E^,  H     be  incident  in  the  half-space  z  >  0  on 
a  perfectly  conducting  screen  in  the  plane  z  =  0,  the  aperture  in  the  screen  being 

denoted  by  S,     Then  the  total  field  in  the  half-space   z  <  0  is 


(8.1) 


X       dz   *     y       az  '        z  dx       dy     * 

H     »  iky  *  |H  ,      H     =  -iku  +  |5  ,   H     ^^  , 
X  ax  *      y  ay  *  *^z      az  * 


where! 

(8.2)  (u,v,w)   -  |-yy(e^,ey,h^)  $  dx'dy' 
and 

(8.3)  xkw--^.^. 

The  functions  e   ,e   ,h     idiich  are  connected  by  the  relation 
X     y     z 

de         ae 
(8.U)  ikh     -  -  — ^  +  ^  , 

z        ax'      ay'  * 

satisfy  the  differential-integral  equations* 

h//\  ^o^'^y    *  ^^#V  lo^'^i^'    =  2nHj(x,y,o), 
(8.5)  ^//^\  lod^'dy'    -  ik^e^  J^dx'dy'    =  2nHJ(x,y,o), 

h/^^x  Jo^'^'    -^  h-^^J  ^o*^'^^'   "  -2nEi(x,y,o), 

when  (x,y,o)  is  a  point  of  S,  If  there  were  no  aperture  in  the  screen,  the  total 
field  would  be  null  in  z  <  0,  but  would  be  E  ,  H  ,  say,  in  z  >  0.  In  the  presence 
of  the  aperture,  the  total,  field  in  z  >  0  is 


(8.6) 


E  =  E°  -  IH  ,  E  =  E°  -  |X  ,  E  =  E°  -.  |H  +  |Z  , 

X  X  az  *  y    y  az  '  m    z  ax  ay  * 

u  tjO  •  i„        Sw     ,,  TjO    .  .  I           aw  „  ,.o        dw 

H_  =    H  -    ikv    -    -r— ,     H  "C    H       +  iku    -    -r-,  H  =    H       -    rr-     , 

X  X  ax'     y  y                    ay'      z  z        az 


+  ^  =  e"^^/R,       R^  -   (x-x')^+(y-y')^  +  z^j  5o  °  ^  "^®"  ^  "  °' 

■«■     Copson  terms  them  "integral  equations"    • 
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My  comment  upon  this  theorem  was  as  fellows'-    -' : 

"There   remains  only  one  question  not  properly  accounted  for  by  this 
analysis,   nattely  whether  or  not  line  charges  along  the  rim  of  the  screen  are 
necessary*.     For  instance  the  proof  of  theorem  h  is  incomplete.     Suppose 
the  equations  (5)   are   solved  rigorously,  under  the  side-condition  (U).     It  is 
not  at  once  evident  -vdiether  the  wave  functions  u,v,w  defined  by  (2)   fulfill  (3). 
Now  it  can  be   shown  that 

Maxwell's  equations  are  satisfied  if  and  omly  if  the  right-hand  integral  van- 
ishes.    Because  e     =  £   ,  e     =  E     in  the  hole,   the  integral  over  S  can  be  trans- 

X    x'  y    y 

formed  into  an  integral  along  the  rim  of  the  hole 

(8.8)  |^/$E^ds, 

vriiere  E     denotes  the  electric  field  tangential  to  the  rim  of  the  hole.     Therefore 
s 

the  condition  for  Maxwell's  equations  to  be   satisfied  is 

(8.9)  Eg  =  0. 

This  is  an  extra  condition  with  regard  to  the  solution  e  ,  e  ,  h  .  For  a  circular 

X       y       z 

hole   the   condition  is  equivalent  tox'e     -y'e     =Oat  the  rim.     This  condition  is 

y  X 

satisfied  in  Copson's  (and  Bethe's)  theory  for  the   small  circular  hole,  and  Copson 
states  that  his  approximate  solution  does  not  violate   (3)»" 

The  validity  of  Copson's  assertions  and  his  theorem  has  not  been  questioned} 
however,   I  claim  that  a  necessary  condition  for  Conson's  theorem  U  to  be   self- 
condistent  is 

(8.10)  e     =  0     on  the  rim  of  S, 

s 

-» 
where  e     denotes  the  projection  of  the  vector  e   =   (e    ,e    )   at  the  rim  upon  the 

tangent  to  the  rim. 

The  condition  (10)   does  not  imply  that  either  e     or  e     or  both  ai^e  finite 

^  y  -1/2 

on  the  rim.     In  fact,   they  generally  become  infinite  on  the   rim  of  order  D         , 

where  D  denotes  the  distance  to  the  rim.     Thus  (10)   does  not  exclude  that  e     is 

•M-  The  same  question  was  raised  by  Bourgin  in  a  review  of  Bethe's  paper  (Bethe'-     -' 
Bourgin    [53] ) . 
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infinitely  large  on  the  rim,  where  e  means  the  projection  of  e  at  the  rim 

upon  the  normal  of  the  rim. 

Incidentally,  e  and  e  should  be  considered  as  limiting  values  of 

the  corresponding  components  of  e  when  the  field  point  tends  to  the  rim  from 

the  interior  of  the  aperture.  The  vector  e  is  not  defined  outside  S  although, 

following  Copson,  we  may  take  it  to  be  zero  there:  in  this  case  (e  ,e  )=(E  ,E  ) 

X  y    X  y 

over  the  whole  plane  z  =  0.  Of  course,  the  limiting  values  of  e  and  e  are  zero 
when  the  field  point  is  tending  to  the  rim  from  outside  S,  In  this  sense,  (lO) 
states  that  e  is  continuous  on  the  rim.  On  the  other  hand,  e  is  in  general 
discontinuous  there. 

To  prove  that  (lO)  is  a  necessary  condition  it  will  be  assumed  that  S 
is  a  finite  part  of  the  plane  z  =  0,  with  a  simple  closed  boundary  curve  s.  More- 
over, it  will  be  supposed  that  e  ,e  ,  and  h  are  absolutely  integrable  over  S 

X     y  z 

and  have  continuous  first-order  derivatives  with  respect  to  x'   and  y'   in  the 
interior  of  S, 

Let  S     be  any  subdomain  of  S  with  boundary  curve  s     such  that  s  and 
s     have  no  points  in  common.     Then  Stokes's  theorem  may  be  applied  to  the  vector 
"e  in  the  domain  S     when  the  field  point  (x,y,z)   is  outside   S.     Thus 


O   '  •*  O 


e  ds, 


where  the  integrand  of  the  surface  integral  is  the  z  component  of  curl  (  C^e). 

Now  using  (U)  and  the  fact  that  (8/dx',  3/dy' )  $=  -(3/ax,3/ay)  $  ,  we  see 

that  the  integrand  is  equal  to 

-ikh  Q-e  -5—  +e  ■^—  • 
z  ^    y  3x    X  3y 


It  thus  follows  that 
"  3x 


(6,11)      -  l:^//   ey^'dy'  +  ^//^^  dx'dy'-ik^h^|  dx'dy'  - /^  e^ds. 

S  ^_  ^_  ^„ 

o  o  o  o 


The  change  in  the  order  of  integration  and  differentiation  is  justified,  since 
all  integrals  involved  are  uniformly  convergent  with  respect  to  x,y,z  and  ab- 
solutely convergent  with  respect  to  e^,  e  ,h  • 
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The   same  remark  applied  if  in  (ll)  the  limiting  case  is  taken  where  S     becomes 

o 

S,  though  it  must  be  implicitly  understood  that  certain  regularity  conditions 
as  to  the  boundary  curves  s^  and  s  are  satisfied.  If  so,  then  (11)  leads  to 

(8-12)  -  i  *  I  -^'™  •  h^i  %<•= 

and  then  (lO)  follows  from  (12)  and  (3), 

As  emphasized  before,  in  the  proof  of  (10)  it  has  not  been  required  that 

e  and  e  are  finite  on  s.  Yet  an  explicit  example  may  be  useful, 
X      y 

If  S  is  a  circular  aperture  of  radius  a,  the  following  functions  are 

admissible  (so  far  as  Integra bility  is  concerned); 

,  2a^-x'^-2y'^      „  \r2   72   T   „  . 
e  =  A        -^— +  B  y/a  -x'  -y'   +  Cy', 

(8,13)  l/^  -^  -y 

^  y _  Pvi 


e  =  A 

y     {12      ,2   ,2 
/a  -x'  -y' 

where  A,  B,  and  C  are  constants.  It  is  easy  to  verify  that  the  corresponding 

function  h  is  given  by 


+  2C, 


(8.1U)  ikh  =  -(3A+B)  2Ll 

^a  -x'  -y' 

2   2    2 
Clearly  these  functions  are  absolutely  integrable  when  x'  +y'  <  a  ,  and  are 

2   2    2 
continuously  differentiable  when  x'  +y'  <  a  .  Unless  A  =  0,  they  are  all 

singular  on  the  rim.  Notwithstanding  this,  e  is  finite  on  the  rim,  as  is 

readily  verified  if  polar  coordinates  are  introduced.  Letting  x'  =  r'  cos  ^f', 

y'  =  r'  sin  ^'  we  have 

'         2   2      ] 

e^  =  JA   Jg-Il!-  *   B\/a^-r'^  [  cos  {?' 

(8.15) 

e  =  -Cr'  -  (2A+B)  \/a-r^  sin  gf' , 

and  so  in  this  case  e  =  -Ca  on  the  rim.  Consequently,  the  condition  (10) 

would  rule   out  the  case  C  /  0,  but  it  does  not  require  A  =  0, 


-  6U  - 

In  view  of  this,  I  shall  now  quote*  from  Copson's  answer:  L^  J 

"Bouwkamp's  argument  depends  on  differentiation  under  the  sign  of  integration 

and  the  use  of  Green's  theorem.     If  e     is  continuous  it  vanishes  on  the  ririi, 

v^ereas,  if  it  is  discontinuous,  it  becomes  infinite  on  the  rim,  as  Bouwkamp 

himself  has  point  out.  And  similarly  for  e  .  Thus  either  e  and  e  vanish 

on  the  rim,  and  (3)  is  a  consequence  of  (2)  and  (U)}  or  else  at  least  one  of 

the  functions  e  and  e  is  infinite  on  the  rim,  and  Bouwkamp's  argument  fail^'« 
X     y 

As  a  matter  of  fact,  my  proof  of  (10)  shows  that  Copson's  assertion  is 

simply  incorrect,  as  is  also  borne  out  by  the  explicit  example  (13) »  Vlhat  matters 

is  that  e  should  not  be  singular  on  the  rimj  this  constitutes  an  intermediate 
s 

case  which  has  been  overlooked  by  Copson. 

Consequently,  I  fully  maintain  my  "criticism"  of  Copson's  paper.  That  is 

to  say,  if  we  intend  to  solve  Copson's  differential-integral  equations  (U),  (5), 

we  should  look  for  those  solutions  e  ,  e  which  satisfy  the  auxiliary  condition 

X       y 

(10)  because  then  and  only  then  will  the  electromagnetic  field  specified  by  (1) 
and  (6)   solve  Maxwell's  equations.     This  information  concerning  the  solution  of  the 
differential-integral  equations  (5),  which  was  omitted  in  Copson's  paper  and 
which  I  therefore   added  in  my  review,   does  not  weaken  the  value  of  Copson's 
theorem.     On  the  contrary,  it  was  meant  to  be  and  is  in  fact  a  further  step 
towards  the  practical  application  of  the  theorem,  especially  in  the  construction 
of  approximate  solutions. 

In  addition,  it  is  now  evident  that  in  a  rigorous  formulation  of  plane 
diffraction  problems  there  is  never  need  of  additional  line  integrals  along 
the  rim.     This  settles  an  old  question  concerning  the  "rigorou^'  extension  of 
the  Huygenft— Kirchhoff  principle  to  electromagnetic  diffraction  problems.     Where- 
as fictitious  boundary  values  of  the  field  vectors  on  an  open  surface  in  general 
require  these  line  integrals  in  order  that  Maxwell's  equations  be  satisfied, 
the  correct  boundary  values  automatically  make  these  equations  vanish  identically. 

At  the  time  of  writing  my  review  of  Copson's  paper,  his  approximate   solutions 
for  the   small  circular  disk  and  aperture  seemed  to  be  correct,  since  the  condi- 


*  Actually,  Copson  discusses  the  complementary  problem,  his  theorem  5.     I  pre- 
fer to  keep  the  argument  in  accordance  with  what  I  criticized  in  Copson's 
work. 
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tion  (10)  happened  to  be  fulfilled  and  Bethe's  earlier  results  were  confirmed. 
However,  mainly  because  of  Meixner's  investigations,   I  have   since  come  to  the 
conclusion  that  Bethe's  as  well  as  Copson's  first-order  solutions  are  correct 
in  the  wave   zone,  but  that  they  fail  in  and  near  the   aperture  or  disk. 

The  fact  that  Copson's  approximate  solution  for  the  small  circular  aperture 
is  incorrect  will  be   shown  by  the   simplest  possible  example,  namely  the  diffi'ac- 
tion  of  a  plane-polarized  wave  impinging  in  the  normal  direction.     In  this  case 
Copson's  equations     (6.U)   -  (6.6)'-     -'  reduce   to 


Uik  ^j:;;^    .     „     „  n.      V,     =  _  li  tL 


(8.16)  e^  =  ^  ^a  -r'^   ;     e^  =  Oj     h     =  -  ^, 


X        Ti     '  '       y         '       2 


\^a  -r' 


if  it  is  assumed  that  the   incident  wave  is  polarized  parallel  to  the  x  axis 

and  of  unit  amplitude   (E^  =  1,  H^  =  1).     In  the   limiting  case  ka  -^  0,  Copson's 

X  y 

equations   (5)   then  reduce  to 

dx/^     z       p  JA     Y         9 


ay 

3x 


where 


2        /         ,^2        /        ,^2        2     2^2 
p     =   (x-x')     +   (y-y')    ,   X  .+y     <  a   , 


If  the  expressions  (16)    are    substituted  in  (17),   the   right-hand  members  become, 
in  the   same  order, 

Oj   -2n  +  nk^(2a^-x^-y^)   =  -2n+0(k^a^)j  -2Tiikx  =  0(ka)j 
thus,   as  in  Copscn's  paper,    (16)   is  an  approximate   solution  of  (17) • 

However,   a  second  solution  is  provided  by 
/u   iQx  Uik  2a^-x'^-2y'^  Uik       x'y'      .   ,  U         y' 

which,  if  substituted  in  (17),  will  make  the  right-hand  members  equal  to 
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-  I  nk^xy  =  O(k^a^);   -2n  +  ^  Ttk^(20a^-3x^-9y^)=  -2n+0(k^a^);  0. 

This  clearly  demonstrates  that  equations  (U)   and  (5)   of  Copson' s  theory 
must  be  handled  very  carefully  in  order  to  determine  an  approximate  solution 
e   ,  e   ,  h   .     In  addition  it  is  to  be  noted  that  the  approximate  solution  (l8) 
is  better  than  Copson' s  (16)   since  in  the  first  case  equation  (17)   is  satisfied 
up  to  terms  of  relative  order  (ka)     and  in  the   second  case  only  up  to  terms 
of  relative  order  ka« 

The  condition  (10)  is  not  decisive  as  to  the   question  vAiether  (16)  or 
(18)   is  the  physical  solution,  since  both  (16)   and  (l8)   are  consistent  with 
(10),     Thus  (10)   is  not  a  sufficient  condition,   at  least  not  for  the  purpose 
of  findino:  the  long-wave  approximation. 

The   solution  of  this   difficulty  is  simple.     A  detailed  investigation  of 

the  electromagnetic  field  calculated  on  the  basis  of  Copson's  (that  is,  Bethe's) 

approximate  solution  (16)   has  revealed  ( Bouwkamp  ••     -^ )  that  the  corresponding 

electric  field,   which  is  throughout  of  order  ka  compared  to  the  magnetic  field, 

is  discontinuous  in  the  aperture.     There  is  no  sense  in  retaining  E     in  the  aperture 

and  ignoring  E   ,   as  Copson  did  in  his  approximation,   since  these  quantities  are 
z 

of  the   same  order  of  magnitude. 

As  was  shown  elsev^ere    ( Bouwkamp ■-     -^ )   the  approximate  solution  (18)   does 
not  lead  to  a  discontinuity  in  the  electric  field  in  the   aperture. 
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IX.  Diffraction  by  Narrow  Slits 

In  this  section  we   shall  be  concerned  with  approximate  solutions  for 
the  diffraction  of  plane  waves  by  narrow  slits.     The  x-axis  is  assumed  to 
be  parallel  to  the  edges,   and  we   shall  consider  only  problems  that  are  in- 
dependent of  X.     That  is,  the  direction  of  incidence  of  the  plane  wave  is 
in  a  plane  normal  to  the  edge.     The  integral -equation  formulation  of  these 
two-dimensional  problems  follows  at  once  from  the  corresponding  formulation 

in  three  dimensions^  if  we  integrate  with  respect  to  the  x-coordinate .     If 

2  2  2  2 

in  three  dimensions  we  have  r  =  (x-x')     +   (y-y')     +  (z-z')     and  in  two 

2  2  2 

dimensions  p     =  (y-y')     +  (z-z*)    ,   then 

00     ikr  /, , 

(9.1)  /     ^    dx'   =  ni  H^^^kp). 

-00 

The  two  principal  boundary-value  problems  for  the  slit  are  to  be  for- 
mulated as  follows: 

Problem  I,  ^  =0  on  the   screen: 

9  =  ^o^'^^   -  ^o^^'-^^   ""  ^i(y'-2)  (z  <  0) 

9  =  ?-L(y,z)  (2  >  0). 

Problem  II,   3^/3n  =  0  on  the   screen: 

^  "  ^o(y*z)  +  9o(y,-z)  -  <i^2^yy-z)  (z  <  0) 

where  the  wave  functions  (ji,   and  ^^  (defined  for  z  >  0  only)   can  be  represented 
in  the   form  of  integrals  extended  over  the  aperture,   the  integrands  contain- 
ing the   aperture  values  of  ^^    and  B^'p/Bn  respectively.     The  two-dimensional 
analogue  of  Rayleigh's  forraulss  are 

t    See  chapter  III. 
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(9.2)      f    {2^(7, z)   =-|^/    ^i(y,0)H^^^ky(y-y)2+z^   )    dy' 

-a 


-a 


vhere  2a  is  the  width  of  the  slit. 

As  in  the  case  of  the  circular  aperture,  we  obtain  differential-integral 

relations  for  9'-i    and  ^^o  ^y  requiring  that  d^/3z  =  3^  /dz  in  the  aperture  in 

the  case  of  proMem  I,   and  $2  =  ^     in  the  aperture   for  problem  II.     Of  course, 

9     is  the  incident  field.     These  differential -integral  relations  are 

(9.3)     -2i  [^  ?Jy,z)4„„  ■=  (k^-  ^)  /  Pi(y',o)H^^^kly-y  |)dy' 

dy       -a- 

a    ra9'2(y'»z')  1  (-,) 

(9.U)        2i  9(^(y,0)   =/  J  -^, V        H^^^kly-y'Ddy', 

-a   j  Jz=o 

where  -a  <  y  <  a. 

Henceforth  we  consider  only  the  case  of  normal  incidence;  that  is,  we 

choose  5?  =  exp(ikz).  We  also  introduce  the  notations 
7  =  8  sine,  y'  =  a8in©',  >:8  =  e 

r.  (e)  =  -  -~-   cos  Q  <fi^  (a  sin  ©,  O) 

F-(«)  =  a  cos  e  I^  ^.(a  sin  e,z)l   . 

Equations   (3)    and  (it)   then  reduce  to 
p  _  /p 

(9.5)  2i  =   (e^  +  a^^  )/        F,  (©' )H^"'-\e  |sin&-sinO' | )   d©' 

dy^   -Ti/2 

"/2       .^^ 

(9.6)  2i  =/    F  (©•)H^^(6|sin9-sin©' |)dO'. 

-Ti/2       ° 

We  now  try  to  find  approximate  solutions  of  the  equations  (5)  and  (6) 
for  the  case  where  e  is  small.  To  this  end  we  recall  that  we  know  the  eigen- 
functions  of  the  corresponding  static  potential  problems,  Kore  explicitly,  the 
eigenfunctions  f(&)  of  the  homogeneous  integral  equation 
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n/2 
(9.7)      /        logV   2|sine-sine' I  WOD  d©'   -  Xf(0) 
-n/2 


J 


are  sin(2n+l)©  and  cos2nO,  where  n  =  0,1,2,,,..    The  relevant  eigenvalues  X 
are  equal  to  -n/m,  where  in  ■  2n,   2n+l  (m  f  o)  respectively;  they  are  zero  if 
m  =  o.     Thus 

(9.8)     -  I  log    |2(sine-sin©')|   -  H  ^  cos2necos2n©' 

00 


*y_    ^^  sin(2n+l)9sin(2n+l)©«. 


[55] 


Both  ©  and  0'  are  between  -n/2  and  n/2.  Equation  (8)  follows  from  the  known 

expansion 

1  °°  1 

-  •^  log|2(cos9'-cos9'' )  I   -  ^  -  cosa^cosmy 

;rfien  we  make  a  suitable  change  of  variables. 

We  need  an  appropriate  expansion  of  the  Hankel-f unction  kernel  in  o\u* 
integral  equations.     This  is  a  series  expansion  in  powers  of  e  in  which  the 
coefficients  depend  on  ©,  6'   and  on  log  e.     This  expansion  is 

H^^^6|sin©-sin©'|)  »  ^  [l^+T^e^+f^^^*  "•  ""V^"  *  O^^^"""^  ^°g  ^^  ]  » 
where 

Yp  =  p  +  log|2(sin©-sinO') I 
(n>o)         ^ 

Y  =  -^^^ — TT   ( sin©- sin©  i)^"  fp  -  1  _  i  ....  i  +  log  |2(sin©-sinO' )  I  ] 
n   „2n/  ,\2  ^  2  n     ^'  '-■ 

2  (nl) 

and  p  is  an  abbreviation  for 

p  =  Y  +  log(  T-  £)-■«•  i    (y  "  0.577  •••  =  Euler's  constant). 

I«t  us  first  consider  problem  II,  Assume  that 
F2(©)  =  fo  +  \^   +  ^2^  *    "'      • 
Then 
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F2(e')H^^^6|8ine-sinO'|)  «  ^   [^0  *  *^^Vl  *  Vo^ 

Substituting  this  in  equation  (6)  gives  us  an  infinite  set  of  integral  equations 

of  the  potential  type,  viz., 

* /f  I  dO'     -  It 
-*      00 

(9.9)^     /f^YpdO.    =-/foV®' 

/f^T^d©'    .  -/fj^d©'    -/f^^T^dO',  etc. 

Since  we  know  from  (7)  or  (8)   that 
/T^d©'    =  n  p, 

we  immediately  get  the   zero-order  approximation 
(9.10)  f^  .  1  . 

If  we  insert  this  result  in  the  right-hand  side  of  the  second  equation  of 
(9)  we  get  an  integral  equation  for  f, ,  namely 


/h\^^'  -SJi-H^^*?)  °°^2©  I 


To  solve  for  f.  Me   assume  f,  has  the  form  f,=  c,+CpCos2e',  where  c,  and  Cp  are 
independent  of  ©',  Then 

y  f-j^f^d©'  =  c^T^dO"  +  c^l^cQs2Q^dQ^   =  c^np  -  •^  C2COS2©. 

Therefore  we  find  by  comparing  coefficients    that 

(9.11)        ^1  '  F"  >  ^  ■"  ^^P""^)  ^°^  ^®  f  • 

>  -* 

On  inserting  f  and  f,  from  (10)  and  (11)  in  the  third  equation  of  (9) 
o     1 

and  evaluating  the  necessary  integrals,  we  find 
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/Vc^®'    =  -  ?/ V^'    -  n?/ V®'-  ^l  '  y  /Y^cos2e'da. 


II  3        1         /■12\^^ 


1(1  ^  ^)cosU©    V  . 


This  in  turn  leads, in  the  same  way  as  before  for  the  function  f, ,   to  an  explicit 
solution  for  f„,  viz.. 

Equivalent  representations  of  the  results  (10),   (11),    (12)   are 


0       p 


(9.13) 


1        1       /I       1  s         2^ 

1  °  H?  -  H  "   ^2  "  II?)   °°^  ® 


f-,    = 

^P 

1 

•2  '  12H 


REMARK.     In  deriving  this  result  we  have  used  various  transformations  of  an 
essentially  elementary  character.     For  convenience  we  list  them  here: 

[sin0-sin6']      =  1  -  |  cos2©-  ^  cos2©'    -2  sinQsinQ' 

[sin»-sin&']^  =  j:  -  '^  cos2e  +  -n  cosU©  +   sin©' (sin3©-6sin©) 

+  sin©sin3©'   +  cos2©' (^cos2©-2)   +  ^osU©' 
/T^d©'   =^     {  -  Up  +   (2p+l)   cos  2©} 

fl  cos2n©'d©'    =  -  #-  cos2n©        (n  >  o) 

/^2'^®'    =  &I  ]  f  P  ■  I  ^   ^I  ■  2p)cos2©  +  |(p  +  ^)cosU©    >► 
y^T^cos26'd©'    "Jslp+r-^  cos2©  -  ^  cosU© 


Finally,   our  result   (13)   can  ice  brought  into  the  following  form:      If  we 
assume  the  field  expansion  in  the   aperture   to  be 
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{9,1k)  3I7  ^'^(y*^)  - 1  LpJ^  -  yW)"-'/' 


then 


(9.15) 


r> 


h       2 

correct  up  to  and  including  terms  of  order  (ka)  /tlogka)  .  The  error,  represented 

by  ...,  is  0[(ka)^. 

It  is  to  be  noted  that  oxir  result  confirms  that  of  Somraerfeld   [^6]  who 
only  gave  terms  up  to  but  excluding  (ka)    ,     A  small  error  of  Sommerfeld  must 
be  noted,  however:     his  A'    should  be  replaced  everyvAere  by  2AV'ii  (see  also 
H6nl  and  3roschwitz  *-     -• ) . 

We  now  turn  to  calculating  the  wave  field  at  large  distances  behind  the 
slit.     Let  p,  6  be  polar  coordinates  with  respect  to  the  positive  z-axis.     Then 


-a 
Inserting  (lU)   and  (15)  we  get 

^2-ly^    ei^^P-3"A)    [l.J(ka)2cos2. 

+  i-(ka)^  J-pcos2e+sin^0-3sin^»+  ^  +  i   [•  +   •••  ]j 


hence 


W2'. 


IpI 


-Rcos2Q+2cos^e+sin^e-3sin^©+^       ^  ^ 


}^' 


where  e  =  ka  and  R  =  real  part  of  p. 

Averaging  over  all  angles,   and  using 


cos  Q     =     sin  ©     = 


cos^     =     sir?0     =  I  , 
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ve  find  that  the  transmission  coefficient  is 

2 


Uelpl 


'¥} 


or,  explicitly, 


(9.16)      r,^-^^-^^ 2  &  *  J(^^^'  ^  2I5  (^-^'*  1 1  * 

rvM-102(ika)1^+  ^  L 


[T.log(ika)]^.  :^ 


Y+log(7-ka) 


+  u      '--^^u: 


3  1         2     n^ 

[Y+log(j4ca)J   +  -jj- 


•••] 


Note  that  this  Tp  is  also  the  transmission  coefficient  for  the  electromagnetic 
case  if  the  incident  wave  is  polarized  perpendicular  to  the  edge  and  the  screen 
is  a  perfect  conductor. 

Problem  I  can  be  solved  in  a  similar  way.  We  assume  F-.(Q')  to  be  of  the 


form 


F, (©')  =  c  cos  0'  +  Cpcos  ©'  +  c^cos  9'  +  •••  , 


where 

c-L  =  gj^  +  gp  e^  +  g^  e^  +  •  •  • 

and  where  g  depend  on  e  (not  on  0')>  but  only  as  some  power  of  log  e.  Then 
^  F^(©')H^-'-^(6|sin©-sin©' 1)  =  giT^cos^O' 

0  0  li  0 

+   e  [gpl  cos  9'+g  ¥  cos  ©'+g,T,cos  ©'  J 
+  6^  [g^T^cos  ©'  +  grY^cos^©'  +  g^Y^cos  ©•  +  g2Y-j_cos  ©« 
+  g^T-LCOs*^©'  *  g^l^cos^Q^   3  +  •••. 

Ve  thus  require  integrals  of  the  form 
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fl  cos^^'d©'. 
-'      n 

The  following  explicit  results  have  been  foiuid: 

yT^cos^O'dO'   =  r  (2p-cos2©) 

yT^cos^e'd©'   =  -g  (3p-2cos2e-  icos^) 

y  Y^cos^O'd©'   -  ^(lop  -  i^062©-  |cosUO-  gcosdo) 

/Y^cos^e'd©'   =  ^2  •{  I  -  3p  +(2p+  i)cos2©  -  ~  cosU©   > 

/Y^cos^e'd©'   =  32   1  I  "  ^^  ""^i^  *  ^)cos2©  -  i  cosU©  -  g^  cos6© 
.     /Y2cos^©'d0'   =  ^  Ap  -  I  +   (||  -  5p)cos2©+(^  ^)cosU©-j^os6©   L 

In  the  process  of  evaluation  all  functions  Y  cos     ©'   were  first  transfonned  into 

a  form  containing  the  cosines  of  the  even  multiples  of  ©  and  ©'   and  the  sines  of 

the  odd  multiples  of  ©  and  ©'    in  a  symmetrical  way,  and  then  (7)   and  (8)  were 

applied, 

2 

To  cope  with  the  operator  6  +  a  — ^  ,  we  note  that 

dy 

a  — ^  cos2©  =  -U 
dy 

a  ^  cosU©  =  l6(2-3cos2©) 
dy"^ 

af^   °   cos6©  =  -12(9-l6cos2©+10cosU©)} 

dy'^ 

thus  we  have 

2 

a^  ^  Ay  cos^e'dO'  =  n 

dy 

2 

a^  1^  /  Y  cos'^©'d©'  =  ^  cos2© 
^,2  >   o  2 

dy 

2 

a^  l-^/Y^cos^©'d&'  =  -^  (2cos2©+cosU©) 

dy 
2 

a^  ^ /y  cos^e'd©'  =  ^  (-2P-1+COS2©) 
dy2-'   1  8 
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2 
&    —  / l^cos^Q' dQ'   =  ^(-6p-3+Ucos2e+  |  cosU©) 
dy 

2d2 
a    — ^ 


/T^cos^e'd©'   -  jfjf??-  -^  -(6p+  |)  cos2e  +  ^cosae} 


Consequently 

2 

2i(e^  +  a^  ^)  /   F^(0' )H^"'"^elsirt9-sin©' |  )de' 
dy 

2 
=  g-L   [l  +  ^(2p-cos2©)  ] 

+  gj    [e^  +  J  e^(2p-cos2©)  J 


+  g,     [|  6^cos2©  +  j6^(3p-2cos2e-  icosU©)  ] 

+  g^   [  ^  6^(-2p-l+co62©)  +  i^  6^  ^  I  -  3p+(2p+  i)cos2©  -  j^osU©  ^  ] 


+  g^    [  I  6^C0S2©  J 

+  g^    [f  6'^(2cos©+cosUe)  ] 

+  «2    I^F  e^(-2p-l+cos2©)  ] 

+  g.     [  ^  6^(-6p-3+Ucos2©  +  ^  cosU©)  ] 

'*  ^1   I^  iIh  ^^  yP"  I  "  ^^P"  |^cos2e  +  icosU©    l] 

+  terms  of  higher  order. 
According  to  the   integral-differential  equation  (5)   this  should  be  equal  to  2i  plus 
terms  of  order  e   •     Therefore  we  find 


^1-1  E^  =  Ys  ^p^-  h'l2^ 


1         £  _      1 

^2  ~  F  "  U  ^6       320 


^U  '  12  ^5  =  ^ 

We  finally  get  the  following  result:  let 


gc  =i  (l-2p). 


INSTITTJ  L..:CES 

NY. 
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(9.17) 
then 


(9.18) 


52f^(y,0)  =  -ika^  C^(l-yW)"'^^^ 


C^  -  1  -  J(p-  |)(ka)2  *  ^(p2  -  7  p  ^  ^)(ka)'*  + 


C,  = 


^  (ka)^  -^  (p-|)(ka)'*  +  ••* 


RiiMARK.  This  result  is  at  variance  with  that  of  Sommerfeldi"  -J  and  that  of 
Groschwitz  and  H6nl'-  -' .  Sommerfeld  gave  (in  our  notation) 


Cj^  =  1  -  ^  (p  -|)(ka)2  + 


Our  result  (17),   (l8)   in  Groschwitz-Honl  notation  would  be 


—   fl 


2i 


+  e 


3?  is  • 


\  -  \  l°eYe  +  I 


I  log  2  .  ^  .  f 


] 


The  transmission  coefficient  for  problem  I  turns  out  to  be 
(9.19)         Ti  -  §2  ^^^^^  &  *  ^(^^^^  ]l  -  f(r^l°g  J  ka)    I 


+   •  •  • 


] 


Note  that   Xi    applies  to  electromagnetic  waves  if  the  incident  electric  field 
is  parallel  to  the  edge. 
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Section  1 

Diffraction  by  an  Aperture 
in  a  Planar  Screen 

by 
N.  Marcuvitz 
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X.     Diffraction  by  an  Aperture  in  a  Planar  Screen  by  Nathan  Marcuvitz 

The  circular  aperture  problem  solved  by  Bouwkamp  can  be   treated  by  a 
somewhat  different,  but  equivalent,  method  that  emphasizes  the  vector  aspect  of 
the   diffraction  problem.     A  dyadic  Green's  function  formalism  is  employed  to 
obtain  a  transverse     vector  integro-differential  equation  for  the  transverse 
electric  field  in  the  aperture.     The  knowledge  of  the  aperture  field  perrdts 
a  direct  calculation  of  the  electromagnetic  fields  at  any  point.     The  desired 
aperture  field  is  obtained  by  solving  first  an  inhoroogeneous  transverse  vector 
partial  differential  equation  and  then  an  inhomogeneous  vector  integral  equa- 
tion of  the   first  kind. 

The  formulation  of  the  problem  is  independent  both  of  the  nature  of 
the  sources  and  the  shape  of  the  aperture.  The  solution  of  the  inhomogeneous 
steady  state  field  equations 

V  X  E  »  ikH  -  M 
^^^••'•^  7  X  H  -  -ikE 

is  to  be  obtained  subject  to  boundary  conditions: 

a)  n  X  E  ■  0  on  the  screen 

(10.2)  ^)Sscat.  '     Hscat,^°     as  r  ->co  (Im  k  -  o>  y;!?  >  O) 

c)  E         ^i/8»  E  ,   ~  —    in  the  s^erture, 

-tan    :''    '     normal         /^ 

where  s  is  the  normal  distance  to  the  rim  of  the  aperture.  As  Bouwkamp  first 
emphasized,  conditions  a)  and  b)  alone  are  not  sufficient  to  insure  a  unique 
field  solution.  Condition  c)  defines  the  singularity  of  the  electric  field  at 
the  rim  of  the  aperture  and,  together  with  a)  and  b),  does  insure  such  a  solu- 
tion. Condition  b)  refers  to  that  portion  of  the  far  fields  due  to  current 
sources  M(r)  (real  or  induced)  at  a  finite  distance  r.  As  Bouwkamp  discussed 
above,  it  is  convenient  to  divide  the  field  region  into  two  parts,  z  <  0  and 
z  >  0,  to  the  left  and  right  of  the  screen,  with  sources  prescribed  in  the 
region  z  <  0,  The  total  field  in  the  region  z  <  0  is  equivalent  to  that  pro- 
duced in  a  half-space,  with  a  perfect  conductor  at  z  «  0,  both  by  the  pre- 
scribed sources  and  by  the  "induced"  magnetic  current  source  n  x  £ 


•«■  i,e.,  parallel  to  the  plane  of  the  screen. 

^^>   M,K.S.  system  out  normalized  so  that  intrinsic  impedance  i/ti/e  of  vacuum  is 
unity. 
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in  the  aperture   (n  «=  unit  normal  vector  in  positive  z  -  direction).     Similarly 
the  fields  in  z  >  0  may  be  regarded  as  produced  solely  by  an  induced  magnetic 
current  density  -n  x  E  in  the  aperture. 

It  is  convenient  to  introduce  a  "half -space"   dyadic  Green's  function 
y(r,r')  defined  by  (cf.[Uo]) 

(10.3)  V  X  (V  X  Y)  -  k^T  -  -ike6(r-r') 

subject  to  boundary  conditions 

a)  Y  ->  0  as    |r-r' |  ^»   oo     (Im  k  >  0) 

b)  n  X  (y  X  Y)  -  0     at  z  -  0, 

where   e  is  the  unit  dyadic  defined  by  e  .  A  -  A,  and  6(r-r')  is  the  three- 
dimensional  delta  function.     Physically  -Y'r,r')    •  e  is  the  magnetic  field 
produced  at  r  by  a  delta-function  magnetic  current  flowing  in  the  direction 
e  at  r'  in  the  half-space.     The  half-space  Green's  function  can  be  defined 
in  terms  of  a  "free-space"   ctyadic  Green's  function  Y^(r,r')  which  obeys 
Eq.   (3)  with  the  omission  of  condition  b).     The  free-space  Green's  function 
is  given  by      (cf.    |U0|) 

(lO.U)  Y^(£»r)  -  -ik(e  +  ^)  g(r,r') 

where   g(r,r' )  is  the  scalar  Green's  function  defined  by 


(10.5) 
Although 


(V^  +  k^)g  «  -6(r-r') 

g  -^  0  as    |r-r' |  ^  oo   (Im  k  >  O). 

ik|r-r'| 


(10.6)  g  -  ® 

Un|r-r' | 

is  a  simple  closed  form,  a  more  convenient  representation  for  subsequent  appli- 
cations will  be  considered  below.     The  half-space  dyadic  Green's  function  Y(r.r') 
is  obtained  by  additive  or  subtractive   superposition  of  two  free-space  dyadic 
Green's  functions  Y^(r,r'),  one  corresponding  to  the   source  at  r'   and  the 
other  to  its  image   at  r'   -  2n  n  .  r'.     In  particular,   for  a  transverse  magnetic 
source  on  the  z  «  0  plane, 

(10.7)  Y(r,r')   =  -2ik(6  +  ^)  g(r,r'). 
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where  r' =  (x',y' ,0). 

The  magnetic  field  due  to  both  the  prescribed  and  induced  sources 
can  be  expressed  in  terms  of  Y(r,r')>  the  field  of  a  "point"   source.     Since 
the  prescribed  sources  are  assumed  to  lie  in  z  <  0,   the  magnetic  field  in 
z  <  0  can  be  represented  by  means  of  superposition  (or  equivalently  by  use 
of  the  vector  Green's  theorem)   as 

(10.8)  H(r)   =  H^(r)   - Jf^^  Y(r,r')-  n  x  E(r«)clS',     z  <  0 

where  the  first  terra  H     is  the  magnetic  field  produced  by  the  prescribed 
sources  in  the   absence  of  the  aperture,   and  the  second  term  is  the  field 
produced  by  the  induced  sources  n  x  E  in  the   aperture.     Similarly,  the 
magnetic  field  in  z  >  0  produced  by  the  induced  sources  -n  x  E^  in  the 
aperture  is 

(10.9)  H(r)   -  +  ff     Y(r,r')    •  n  x  E(r')dS',     z  >  0. 

It  is  evident  from  Eq.  (2)  that  the  field  representations  (8)  and  (9)  satisfy 
the  field  equations  (1)  in  z  <  0  and  the  boundary  conditions  (2a)  and  (2b), 
Moreover,  since  n  x  (V  x  Y)  has  a  jump  discontinuity  of  value  -ike.6(p-p') 
at  z  "  z»  "  0,  it  follows  that  Eqs.  (8)  and  (9)  yield  values  of  n  x  E(r) 
that  are  continuous  as  z  ->  ±  0  and  equal  to  the  value  of  n  x  E  in  the 
aperture.   The  requirement  of  continuity  of  the  H(r)  x  n  given  by  (8) 
and  (9)  in  the  aperture  repion  (z  =  *  O)  imposes  the  condition 

(10.10)  -H  (r)   x  n  =  2/T^  n  x  Y(r,r')   x  n  .  E,  (r')dS',  r  ^  (x,y,±0) 

-0  -         ^      J  J  ap  -  -  -  —      ~t,    - 

"in  the  aperture" 

on  the  transverse  electric  field  E.    in  the  aperture.     The  continuity  of 
nxYxnatz=0  should  be  noted;   however  z  is  not  permitted  to  equal 
z'   =  0  since  the  integral  in  (lO)  becomes  divergent  at  r  =  r ' ,     Because 
of  this  latter  fact,   as  Bouwkamp  has  emphasized,     Eq,   (10)   is  not  a  true 
integral  equation  but  may  be  called  a  pseudo-integral  equation  for  E, (r). 
In  view  of  the  representation  (7)>  Eq.   (10)  may  be  rewritten  as 

nxV.V.xn         /y 

(10.11)  -H  (r)  X  n  -  Uik(6,    -         l         )    •//    g(r,r' )E.(r' )dS' ,  r  in  the  aperture. 


In  a  rectangular  x   ,  y  coordinate   system,   5(p-p')   =  6(x-x' )5(y-y ' ) » 
6.    "XX     +yy     =  transverse  unit  dyadic. 
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where  the  interchange  of  differentiation  and  intep-ration  is  permissible  since 
the  integral  exists  before  and  after  interchange.     Since  g(r,r')  is  integrable 
even  In  the  liirdt  z  »  z'   =  0,  Eq.    (11)  is  a  true  integro-differential  equation. 
To  emphasize  this  fact  we  may  rewrite   (11)  in  the  form*: 

(10.12a)  curl?  F  -  k^F  =  -ikH  (p)  x  n 


(10.12b)  F(p)   =  U   rr     g(p,p')E,(p')dS' 

~    -^  J  J  ap       ^  ~        u  —■ 


where  we  have   defined 

(10.13)  curl?  F  =  n  X  V. V,    x  n   •  F  »  V.y.    .  F  -  V?  F, 

and  where     p         =   (x,y,0)  is  the  coordinate  vector  in  the  aperture.     It  is 
necessary  first  to  obtain  the   general  solution  of  the   vector  differential 
equation  (12a)   for  F,   and  then  to  solve  the  integral  equation  (12b)   for 
E. ,     The  arbitrary  constants  in  the  resulting  solution  for  E.    are   determined 
by  imposing  the  remaining  boundary  condition  (2c). 

A  general  solution  to  Eq.   (12a)   can  be  obtained  in  a  variety  of 
ways,   depending  on  the  nature  of  the  excitation  H     and  the   shape  of  the 
aperture.     It  thus  appears  desirable  to  specialize   at  this  p^nt. 

Diffraction  of  a  Plane  Wave  by  a  Circular  Aperture 

Since   an  arbitrary  source  distribution  can  be  resolved  into  plane 
wave  constituents,  it  is  basic  to  consider  a  plane  wave  incident  on  the  aper- 
ture.    There  are   lwo  independent  types  of  vector  plane  waves:     the  E-  and  the 
H-  waves-;H4^,   distinguished  by  their  polarization.     If  a  rectangular  coordinate 
system  with  origin  at  the  center  of  the   aperture  is  oriented  so  that  the  plane 
wave  is  incident  in  the   xz-plane,    the  transverse   field  distribution  of  an  E-mode 
wave  is  given  by    ^    ■» 

•«•  Bouwkamp  writes  Eq.   (12a)   in  a  somevrhat  different  form.     From  (12a)   one  notes 

that 

(i)  -k^V^    •   F  =-ikV^    •  H     x  n  =  -k^E 

t      ~  t       ~o       ~  on 

whence  on  expansion  of  curl     of  (12a) 

V  y 

(ii)  {t.    +  k^)  F  -  ik(6,    +  -V^)    .  H     X  n  =  |L  E^^    . 

X,  ^  X         .£.  -vO       'v.       dz  ^ot 

Eqs,  (i)  and  (ii)  which  together  are  equivalent  to  Eq.  (12a),  are  employed  by 
Bouwkamp  but  with  ^  replaced  by  F  x  n, 

•M-*  Cf.  U$J ;  Section  26  contains  a  description  of  a  complete  orthogonal  set 
of  vector  plane  waves. 
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ik  X 

X     ik  X 

(lO.lUa)       e'(p)  =  h'(p)  X  n  =  -V^  ^- =  -e  ^  x 

and  that  of  an  H-mode  wave  is  given  by 

ik  X    , , 
X     ik  X 

(lO.liib)      e"(p)  =  h"(p)  X  n  =  -V.  x  n  ^^r—   =  e  ^  y^  , 

where  e  and  h  denote,  respectively,  the  transverse  electric  aiid  magnetic  field 
distributions  of  the  mode  in  question.  With  Bouwkamp's  choice  of  plane  wave 
excitation,  the  unperturbed  transverse  magnetic  field  H  (p)  x  n  in  the  aperture 
is  a  superposition  of  an  E-  and  H-  mode: 

(10.15)  H^(p)  X  n  =  I'h'(p)  X  n  +  PWCp)  x  n, 

where  I'    and  I?'    are   twice  the   incident  amplitudes  of  the   transverse  magnetic 
fields  of  the  plane  waves  defined  in  Eqs,    (lU) .     For  Bouwkamp's  choice  of  plane 
wave  excitation,* 

(10.16)  I'   =  -2cos0  ,  T"   =  2sinQ)  cose   ,  k     =  k  sin©     , 

'^O'  '^O     0'   X  0 

The  solution  of  Eqs.  (12)  for  a  composite  wave  of  the  form  (15)  is  un- 
necessarily complicated  since  the  symmetry  of  the  excitation  is  concealed. 
Accordingly  it  is  desirable  to  represent  F(p)  as  a  superposition  of  an  E-  and 
H-  mode  component,  viz: 

(10.17)  F(p)  =  ri'F'(p)  +  r'F"  (p)  "I  ; 
then  (I2a)  decomposes  into 

(10.17a)      curl^F'  -  k^F'  =  -ikh'(p)  x  n 
(10.17b)      durlfp  -  k^P'  =  -ikh«'(p)  x  n. 


and  correspondingly  (12b)  decomposes  into 
(10.18a)       F'(p)  =  U  /T'  g(p,p')E'(p')dS' 

(10.16b)       P'(p)  =  h  ff   g(p,p')E"(p')dS', 

J  J  ap 
where 

(10.l8c)       E.(p)  =  I'E'(p)  +  r'E"(p). 
^t  .V     -^    ^  -^    ,^ 


♦  See  Bouwkamp's  definition  of  the  angles  ^  ,  ©  ,  p  on  p.itY 
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Solutions  of  the  vector  partial  differential  equations   (17)  can  be 
obtained  readily.     In  fact  the  general  solution  of  Eqs.   (17)  may  be  expressed 
in  the   form: 

(10.19a)  ^fl  =  felf>^  .  r  2%V  e(-^  -i^  "<?) 

k  m=l 


00  J  (kp) 

+  y  2a'V^x  n{—, cos  m0) 

'=—      m  t    "       k  ^ 


m=o 


(10.19b)  ^  =  ^^ie>n  ^  g  2p  V.x  n(  JH^  cos  mj2) 


ik         ,2,2  -^     '^m  t 

k  -k  m=o 

X 


00  J  (kp) 

*  H  ^PJjV^x  n( -HL sin  m9)    . 

in=l 

From  Eqs.   (13)   and  (llO  we   see  that  the  first  term  in  Eq.   (19a  or  b)  is  a  parti- 
cular solution,   while  the  remaining  "H-mode  terms",  which  alone  can  satisfy  the 
homogeneous  equations   (17),  represent  the  complementary  solution.     The  ^-polar- 
ization    of  the  complementary  solution  is  of  the  most  general  form;  however,  it 
may  be  delimited  by  using  the  symmetry  properties  of  the  field.     In  view  of  the 
rotational  symmetry  of  the  structure   about  the  center  of  the  circular  aperture, 
the  ^-dependence  of  the  field  is  determined  by  the  nature  of  the  excitation.     The 
symrietry  of  the  excitation  is  evident  when  a  solution  in  powers  of  ik  is  considered. 
As  Bouwkamp  has  shown,   only  the  E-mode  component  (l9a)  contributes  to   zero-order 
in  ik,  v/hereas  both  the  E-  and  H-  components  of  F  contribute  to  first-order  in  ik. 

To  make  explicit  the  perturbation  solution  in  powers  of  ik  one  employs 
the  scheme: 

(10.20)       h(p)  =  h^(p)  +ik  lL,(p)  +  (ik}\(p)  +  ••• 


F(p)  =  F^(p)  +ik  F,(p)  +  (ik)%(p)i 


then  Eqs.  (17)  decompose  into 

(10.21a)       curl^F  =  0 

Xr'O 

(10.21b)       curl^F,  =  -h  X  n 

t~l  -o  ^ 

(10.21c)       curl^F^  +  F  =  -h,  x  n 

t~2  --o    -"1  ^ 

(I0.21d)       curl^F^  +  F.  =  -h„  x  n 

t-xj  '-I       -^d              'J 


t  Note  that  x  =  p  cos  ^  \     y  =  P  sin  {?. 
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(the  superscripts   «   and  "   are  omitted).     It  should  be  noted  that  the  zero-  and 
first-order  solutions  ^^  ^"^"^  ^1  ^®  ^'^^  determined  solely  by  Eqs,   (21a  and  b); 
to  determine  the  transverse  divergence  of  these  two  solutions,  Eqs,   (21c  and  d) 
respectively  are  required.     Use  of  Eqs,   (lU)  in  Eqs,   (21a  and  c)  indicates  that 
in  zero-order  both  F     and  the  corresponding  aperture  field  E  (p)  have  no  pre- 
ferred dependence,   i.e.,  they  are  radial.     Furthermore,  Eqs.   (lU)   reveal  that 
the  electric  field  (or  b  x  n)  of  an  incident  E-  or  H-  mode  is  polarized  along 
the  X     or  y     directions  respectively.     In  accordance  with  this  symmetry  Eqs, 
(21b  and  d)  indicate  that  in  the   first  order  the  corresponding  aperture  electric 
fields  must  have  the  follovdng  polar  components^: 

(10.22a)  E'^cos  9>  E"  -^  sin  5? 

E '  /v^sin  ^  ^3"^  *^°^  ^  • 

From  (22a)  the  rectangular  components  follow  by  the  transformation 
(10.22b)       E  =  cos  9>  E  -  sin  ^f  E^ 

X  p  {^ 


\ 


E     ■  sin  or  E       +  cos  0  E_,  , 

It  is  convenient  to  deal  with  rectangular  field  components  since  Eqs.   (18)   imply* 
that  the  angular  symmetry  of  the  rectangular  (but  not  polar)  components  of  E  and 
F  are  similar.     Since  the   symmetry  of  the  first-order  aperture  field  is  known, 
we  can  give  a  more  detailed  characterization  of  the  rim  singularity  than  that 
gi^en  in  (2c),     Equations   (22a)   state  that  in  the   first  order  the  rim  condition 
(2c)  may  be  decomposed  into 


(10,23) 


/s  E'  '>-  cos  9*        ^a'^^  ^^'^  ^ 
/s  W'^  sin  (jl         E"'^/s  cos  (jl  » 


Bouwkamp  has  pointed  out  that  the   zero-order  solutions  obtained  in 
the  literature s-»H<-  have  usually  been  correct,  but  that  this  is  not  the  case  for 
the  first-order  solutions.      Since   the  difference  between  the  method  employed  here 
and  that  of  Bouwkamp  is  most  evident  in  the  first  order,  only  this  order  of  solu- 
tion will  be  considered  below, 

T    See  note  on  previous  page, 

•«•     Since  g(p,p')   is  independent  of  the  ^  orientation  of  the  coordinate  axos. 

*♦  See  section  VII.     For  a  derivation  employing  the  methods  herein    ,    cf.      the 
author's  "Coupling  of  Waveguides  by  Small  Apertures",  p. 68,   rej  irt  Rl,.?-^!?, 
PIB-106,  Polytechnic   Institute  of  Brooklyn. 
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First  order  solution  in  (ik). 

The  general  solutions  (19)  permit  a  ready  determination  of  the  first- 
order  form*  of  F(p).  It  is  relevant  therefore  to  consider  at  this  point  the 
first-order  form  of  the  integral  equations  (18),  Let 

g(p,P')  -  g„(p,P')  +  (ik)g  (p,p«)  +  ...  =  — i +  ^  +  ... 

(10. 2U)  Un|p-p'| 

E(p)  =  Ep(p)  *   (ik)E^(p)  +  •••; 

Then  in  view  of  the  corresponding  expansion  of  F(p)  in  (20)   and  the   stated 
radial  nature  of  E^(p),  one  has  the  first-order  integral  equation  (for  either 
the    '  or  "  component): 

where  g^(r,r')   is  the  static  form  of  the  Green's  function  in  (5),   defined  to 
within  a  constant  by 

V^g  (r,r')   =  -S(r-r'). 

In  an  oblate  spheroidal  coordinate  system    6,  ^  associated  with  an  aperture  of 

radius  a, 

X  =  a  cos  6  cos  ^ 

(10.26)  .     ^     .     ^ 

y  =  a  sin  9  sm  ^, 

where  0  <  ^  <  2ti  and  0  <  6  <  ti/2.     The   static  Green's  function  g     can  be  represented 
diagonally  in  this  coordinate  system  in  terms  of  a  complete   set  of  orthogonal 
functions  P  (cos©)   ^^^n  which  are  periodic  in  ^  and  vAiose  derivative  with  respect 

to  0  vanishes  at   6  »  -r.     The  latter  property  laplles  thwt  n  pnd  m  are  either  "both 

even  or  both  odd  integers,  and  that  m  <  n.     The  desired  diagonal  representation 

is 

1      00      n  Fc        M  "l2 

^o^P'?.'^   -k^T.  \(2n*l)  tHT^  ^^°M  P^(c=O30)P^(cose.)cosm(5?-^.) 
7X  n=o  m=o  L  J 


.[66][wl 


(10.27) 

where  e  is  the  Neumann  number  and  equals  1  or  2  depending  on  whether  m  =  0  or 
>  0.  On  substitution  of  (26)  and  (27)  one  obtains  as  the  "diagonal"  form  of  the 
integral  equation  (25): 


vr  The  first-order  solution  can  likewise  be  obtained  by  means  of  Eqs.    (21). 
■iHJ  Note  that  Eq.    (27)   represents  a  convenient  form  for  some  of  Bouwkamp's  integral 
theorems . 
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(10.28)  cosm(9f-;2')>, 

where 

6(e)  =  cos  e  E,(p), 

It  is  evident  that  Eq.  (28)  can  be  solved  immediately  on  representation  of  F-,(p)ln 
a  series  of  ^(^0  5©)^^^^^  functions.  Two  cases  will  be  distinguished,  one  in 

which  F,  =  F •  and  the  other  in  which  F,  =  F" . 

•vX        ^±  -wX        />^l 

E-mode  Solution 

In  view  of  Eqs.   (lUa)   and  (22a)  we  see  that  the  only  terms  of  the 
general  E-mode  solution  (19a)   that  yield  contributions  to  the   first  order  in 
ik  are  2  2 

(10.29)  J   (kp) 

+     2aJJ.T  n(— t: 9in30) 

where  the  omitted  a  ,  a'  terms  have  the  wrong  symmetry  and  the  omitted  ik  xx 

m'     m  h     J  J  X  ~o 

tenn  has  been  considered  in  the   zero-order  approximation.     To  the  desired  order 

Eq.,     (29)  becomes 

2  2  ~1 

(10.30)  Fj(p)   =    L-  i  +(sin^6^-a^k^)^  +   ((a-L+a^)k^+2sin^e^)^os2i?  x^ 

j^2  2 
+   (a^^a^)-^  sin252  y^ 

where  we  have  put  k     =  ksin©   .     Substituting  p  =  a  sin©  and  noting  that 

(10.31)  sin^©  =  I   fl-P^Ccos©)]   =  ip^Ccos©), 

F^Cp)-    \a^-  ^  -(sin^e^-a^k^)| 
(10.32) 


F|(p)=   ra3_-  \  ^(sin^e^-a^^k^)^  +    (a^k^-sin^e^)^P2(cose) 

?  9   9 

2  2       a       2  "1  k  a     2 

+  ((a^+a^)k  +2sin  ©Q)2r  P2(cos©)cos2$^  x  +(a^-a^)-^j— P2(cos©)sin2{?y^  . 


The  integral  equation  for  the  first-order  aperture  field  E'(p)  is  given  by  (2^5) 
with  (32)   as  the  left-hand  member.       On  equating  coefficients  of  corresponding 
terms,  one  finds  that 

(10.'^3)  coB«  E-[(p)   =    [A^+A2P2(cos©)+A2P2(cos©)cos2dx^+B2P2(cose)sin2{2y^, 

where 
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2 

IS. 

2 

I  naA^  =  ((a^+a  )k^  +  2sin^©^)  |j- 
3  „„T32   /     \  k  a 


The  application  of  the  rim  boundary  condition  (23)  to  (33)  yields 


2    2 


A 


'2 

T 


A  -  ^  -3A"  =  Oj 


4 


from  this  the  two  arbitrary  constants  a,  and  a.  follow: 


2       2 
a-,k  =  -  sin  6 
3  o 


,k^  =  1 


(to  first  order). 


Equation  (33)  then  yields  for  the  first-order  E-mode  component  of  the  electric 
field  in  the  aperture: 


(I0.3ii)  E-[(p)  =  j^ 


H-mode   solution 


/a  -p^   (2-sin^©^)x^  +  °     (l+sin^9  ) 

-p 


/a  -f 


As  in  the  case  of  Eq.   (29)   the  only  terms  in  the  general  H-mode 
solution  (19b)  which  are  of  interest  for  the  first-order  solution  are: 


(l-kV/2) 
(10.35)  F.(p)=         /.         y.-2p,V,X  n(^i^  cos0k2e3V  n(l:^  cos 30). 


1,2  ^2 
k  -k 

X 


^Jo^^'^l^'t' 


k 


where   the  coefficients  .8^  and  p^  have  been  employed  to  remove  in  first  order 
the  ikx  y^  term  of  the  particular  solution,   and  where   all  p     and  p •    terms  vdth 
the  wrong  symmetry  (cf.  Eqs.   22)   have  been  omitted.     Using  the  relations   (3I), 
we  can  rewrite  Eq.    (36)   to  the  first  order  as  follows: 

2   2 
J{(p)   =  -(Pi-^P3)^^  P2(cos©)sin252x^  + 


(10.37) 


^  -Pl-(Plk2-tan\)  ^ 

P  P  2  ^2  -\ 

-O^k  -tan  Q^)^  P2(cos&)  +  ((3^-8^)k^-2tan^9^)|j^  P2(cose)cos2i?  y^. 


As  before,  if  we  substitute   (37)    into  the  integral  eqaation   (28)   for  the   aperture 
field  Elj^  and  equate  the  corresponding  coefficients  we  obtain 
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co8«  1^'  (o)  -  C2P2(cose)sin29)x^+   rD^+D2P2(cos©)+D2P2(cosO)cos20  ly^   , 

(10.38) 

2  2 
vhere  3       -2         /„     „   >   k  a 

H  "^^2  =  -(i^i^pj)  -2ir  2 

X 

2 

J  naDg  =  (tan^©^-p^k^)  ^ 

2 
^  naD^  =  ((p^-33)k^-2tan\)  |j^       . 

Furthermore,  the  imposition  of  the  rim  conditions  (23)  on  (38)  leads  to 


°o 

^2 

4- 

2 

3D2  = 

4 

he 

arbitrary  < 

constants 

p^k^ 

=  -tan^e 
0 

2  2 

p,k     =  sec  ©         (to  the  first  order). 

Equation  (38)   then  yields  for  the  H-mode  component  of  aperture  field  in  the 
first  order: 

(10.39)  E"(p)   =  -  4-    12  /a'-p'^  y^    "  ° 


=..  (p)  =  -  1^  [j  v/?^ 


.0         ^2 — 2 


In  view  of  Eqs.   (I6)    and  (I8c)  we   find  by  superposition  that  the 
total  first-order  aperture  field  for  the  Bouwkamp  choice  of  incident  wave  is 

i^l(£^=  IT   [v^^  (2I.x^.I"y^-I.sin\x^) 

(1°-'^^)  ^-^—  P    •    (l'x.I"y     .I-'sin^O 

"^  ^o      ~o  o 


-p 

or,   as  Bouwkamp  obtained. 


x)l 


2 

at    — 

3 
(lO.UOb) 

ft 


L      /a^.p^  (-iic2H  X  11  -  V  E      ) 
n   I  f        "^  —o    -rf         t  on 


+  — p   •    (-ikH  xn  +  V,  E     )    I 

^  -o    -w        t  on      I 


ri    2 

/a  -p 

where  H     and  V.E       are  the  unperturbed  values  evaluated  at  the  center  of  the 
~o  t  on 

aperture . 
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